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TERRAIN  ANALYSIS  BY  COMPUTER 

A.  Keith  Turner  and  Robert  D.  Miles 

Purdue  University 

INTRODUCTION 

Geomorphology  is  becoming  an  increasingly  quantitative  science. 
Many  geomorphologists  now  view  the  landscape  as  a  series  of  open 
physical  systems  tending  toward  energy  equilibrium  (l).   Such  scien- 
tists naturally  attempt  to  describe  terrain  in  unambiguous  meaningful 
quantitative  terms.  Accurate,  reproducible,  numerical  descriptors  of 
terrain  have  become  of  great  concern  to  several  scientific  disciplines; 
including  engineering  geology,  military  geology,  soil  science,  and 
land  locomotion  engineering.   This  paper  describes  techniques  for  the 
analysis  of  terrain  using  the  computer.   Further  additions  and  refine- 
ments to  the  techniques  are  anticipated. 

This  research  was  sponsored  by  the  Indiana  State  Highway  Com- 
mission and  the  Bureau  of  Public  Roads;  however,  this  paper  has  not 
had  the  benefit  of  their  review.   The.  opinions,  findings,  and  con- 
clusions expressed  in  this  publication  are  those  of  the  authors  and 
not  necessarily  those  of  the  Bureau  of  Public  Roads  or  the  Indiana 
State  Highway  Commission. 

The  authors  also  wish  to  acknowledge  the  assistance  of  Dr.  V.  L. 
Anderson,  Department  of  Statistics,  Purdue  University,  who  reviewed 
the  statistical  implications  of  this  study. 


TERRAIN  SAMPLING  PROCEDURES 

Map  digitizations  the  reduction  of  graphical  to  numerical  data, 
is  required  before  computer  analysis  is  possible .  Development  of 
simple,  logical,  and  reproducible  sampling  procedures  is  an  important 
prerequisite  to  satisfactory  terrain  analysis.   Sampling  can  be  per- 
formed at  a  series  of  irregularly  spaced  points  or  at  a  regular 
spacing  forming  a  sample  grid. 

Irregularly  spaced  samples  can  often  describe  the  terrain  better 
with  equal  or  smaller  numbers  of  sample  points  than  can  grid  systems 
which  tend  to  sample  repetitiousiy  in  uniform  areas.  Irregularly 
spaced  samples  are  favored  by  surface  fitting  operations,  such  as 
trend  surface  analysis,  since  the  smaller  number  of  sample  points 
minimizes  computational  complexities. 

With  gridded  samples,  the  position  of  a  sample  in  the  array 
automatically  denotes  its  X  and  Y  coordinates.   This  feature  may  ease 
data  storage  problems  within  the  computer  and  can  be  used  to  advantage 
when  computations  require  knowledge  of  adjacent  points,  since  lengthy 
data  search  routines  are  unnecessary.  While  data  are  most  frequently 
collected  on  rectangular  coordinate  systems ,  polar  coordinates  have 
been  used.  Stone  and  Dugundji  (ll)  discuss  the  merits  of  various  systems, 

This  research  used  gridded  data  to  analyze  the  roughness  and  vari- 
ability of  homogenous  areas.   These  values  in  turn  allowed  the  deter- 
mination of  appropriate  irregularly  spaced  sampling  distributions  of 
larger  areas.   Gridded  samples  were  also  used  to  develop  hypsometric 
integrals  of  various  terrain  types . 


DESCRIPTIVE  MEASURES  OF  HOMOGENOUS  TERRAIN 

Two  techniques  for  computer  analysis  of  homogenous  terrain  types 
have  been  evaluated.  One  of  these  measures  the  surface  roughness;  the 
other  is  a  technique  for  estimating  the  hypsometric  integral  whereby 
the  labor  of  calculating  this  measure  is  substantially  reduced. 

Estimation  of  Surface  Roughness 

A  technique  for  measuring  roughness  has  been  suggested  by 
Hobson  (k) .   This  study  has  used  a  modified  version  of  his  program 
VECTOR . 

Description  of  Program  VECTOR 

A  rectangular  array  of  elevation  readings  is  the  basic  input  data. 
A  set  of  intersecting  triangular  planar  surfaces  are  defined  by  groups 
of  three  adjacent  elevation  readings  (Figure  1A).   Two  different  sets 
of  triangles  can  be  obtained  from  the  same  data  by  redefining  the 
triangle  corners  as  shown  in  Figure  IB.   Type  1  triangles  have 
"northeast -southwest"  diagonals;  while  Type  2  triangles  have  the 
opposite  diagonals. 

Normals  to  these  planes  are  represented  by  unit  vectors.  Mean 
vector  orientation,  vector  strength,  and  vector  dispersion  are  computed 
using  methods  defined  by  Fisher  (3)  and  described  by  Watson  (15). 
Vector  strength,  El,  is  obtained  by  using  the  direction  cosine  method  (5). 
The  standardized  vector  strength,  R,  (where  R  equals  Rl  divided  by  the 
number  of  triangles  J  ranges  in  value  from  zero  for  no  preferred 
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orientation  to  one  for  identical  orientation.   Fisher's  dispersion 
factor,  K,  (3)  indicates  the  variability  or  spread  of  the  unit  vectors; 
it  takes  on  small  values  for  highly  dispersed  distributions  and  extremely 
high  values  for  low  dispersions. 

"Smooth"  areas  generally  have  high  vector  strengths  and  low  dis- 
persions (high  K  values)  as  shown  in  Figure  1C.   Such  areas  may  be 
flat  or  may  have  a  regional  tilt-   "Rough"  areas  of  non-systematic 
elevation  changes  yield  low  vector  strengths  and  high  vector  dispersion, 
thus  low  K  values  (Figure  ID). 

Description  of  Test  Sites 

Twenty- five  sample  areas  were  used;  twelve  of  these  were  located 
in  Indiana.   Eleven  sites  were  measured  on  Is  2^, 000  scale  topographic 
maps;  the  twelfth  was  located  en  a  1;62,500  scale  map  and  corresponds 
to  one  of  the  1:2U,000  scale  test  sites.   These  sites  represent  several 
characteristic  Indiana  landforms  including  karst  plains  (2  sites); 
hills  (3  sites),  plateaus  (l  site),  escarpments  (2  sites),  ground 
moraine  (2  sites),  ridge  moraine  and  dunes  (l  site),  and  outwash  plains 
(l  site).  A  thirteenth  site  was  selected  to  analyze  drumlin  topography; 
part  of  the  Weedsport  Quadrangle  (l;2U,000) ,  New  York  State,  was  used 
for  this  purpose.   Twelve  micro-terrain  maps  were  selected  for  analysis 
from  a  series  of  specially  surveyed  maps  prepared  by  Stone  and  Dugundji 
(ll).  These  maps  were  prepared  with  contour  intervals  ranging  between 
one-tenth  and  one  foot  and  scales  from  1:180  to  1;4800.  A  variety  of 
California  terrain  types  is  represented  including  flood  plains,  badlands, 
and  a  variety  of  desert  landforms  (playas,  pediments,  wadis,  dunes). 


Results  of  Terrain  Analysis  by  Program  VECTOR 

Analysis  of  twenty- five  test  sites,  encompassing  a  variety  of  land- 
forms  and  map  scales,  by  program  '/ECTOR  suggests  that  meaningful  terrain 
analysis  can  be  performed  on  the  computer.   The  following  summarize  the 
results  of  the  current  analyses; 

1)  Program  VECTOR  is  equally  suitable  for  analysis  of  micro  and  macro 
terrains;  provided  that  maps  of  suitable  scale  and  contour  intervals 
are  available  for  the  test  sites . 

2)  Fisher's  dispersion  factor,  K,  is  an  excellent  descriptor  of 
terrain  roughness.  Roughness  is  described  as  "the  presence  of  non- 
systematic  elevation  changes"  (see  Figure  1C  and  ID).   Fisher's 
dispersion  factor  is  defined  as 

K  ■  H  "  1 

N  -  Rl 

where  N  is  the  number  of  observations  (triangles) s  and  Rl  is  the  vector 
strength  (3*k). 

For  smooth  surfaces,  all  the  unit  vectors  tend  to  the  same 
orientation;  thus  Rl  approaches  N,  and  K  -  Rl  tends  to  zero,  and 
accordingly  K  approaches  infinite.   Conversely,  for  very  rough  terrains 
Rl  tends  to  zero,  and  thus  K  approaches  one.   The  smallest  K  value 
obtained  in  this  study  was  31°   This  was  measured  for  a  wadi  in  the 
Mohave  Desert  (ll,  Figure  29)  while  the  largest  K  value,  371+>773»  was 
measured  for  a  playa  surface  (ll,  Figure  55).   Indiana  sites  gave  less 
extreme  values.   The  roughest  Indiana  site  was  part  of  the  New  Albany 
Escarpment  (K  =  8l)  while  the  smoothest  site  was  the  ground  moraine 
just  west  of  Tipton  (K  =  30,  OW).   Because  of  the  very  large  range 
of  K  values,  the  log  K  transformation  is  frequently  useful. 
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3)  All  test  sites  were  analyzed  by  program  VECTOR  using  Type  1  and 
Type  2  triangles  (Figure  IB).  No  important  differences  were  observed 
in  any  of  the  measures  as  a  result  of  type  of  triangles  used  in  the 
analysis.  Figure  2  shows  the  values  of  log  K  for  Type  1  triangles  plot- 
ted versus  log  K  for  Type  2  triangles;  all  points  lie  close  to  the 
k5     line. 

k)     Since  K  is  a  measure  of  vector  dispersion  it  may  be  correlated  with 
measures  of  the  variability  of  attitudes  of  the  triangular  planes.   The 
log  (variance  of  dips)  was  plotted  against  log  K  (Figure  3B)  and  a  least 
squares  linear  regression  line  fitted.   The  linear  correlation  coefficient, 
r  =  ~0.96k,    indicates  a  strong  correlation  between  log  (variance  of  dips) 
and  log  K. 

If  log  (maximum  dip)  is  plotted  against  log  K  (Figure  3A)  and  a 
least  squares  regression  line  fitted,  r  =  -0.950.   The  difference  between 
the  two  correlation  coefficients  is  not  statistically  significant.  A 
strong  correlation  between  log  (maximum  dip)  and  log  K  is  geomorphically 
reasonable;  since  steep  slopes  are  normally  associated  with  rough  topog- 
raphy.  However,  given  two  areas  with  the  same  maximum  slope  value,  one 
might  have  a  much  greater  range  in  slope  values,  and  hence  a  larger 
variance  of  dip  and  greater  roughness.   Thus  maximum  dip  may  not  be  quite 
as  good  a  descriptor  of  terrain  roughness  as  variance  of  dip. 

Initially  it  was  believed  that  a  similar  correlation  might  exist 
between  log  (variance  of  strikes)  and  log  K.   Figure  3C  is  a  plot  of 
log  (variance  of  strikes)  versus  log  K.   The  linear  correlation  coeffi- 
cient r  =  0.009  indicates  essentially  no  linear  correlation  exists 
between  log  (variance  of  strikes)  and  log  K.   This  may  be  due  in  part  to 
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the  test  sites  all  having  highly  variable  strike  values .   Careful  exam- 
ination of  the  results  suggests  that  K  may  be  related  to  a  function  of 
the  variabilities  of  both  strikes  and  dips.  At  the  present  time, 
however,  a  satisfactory  model  has  not  been  discovered,  although  several 
have  been  tried. 

5)  The  orientations  of  the  triangular  planar  elements  become  progres- 
sively poorer  estimators  of  the  true  terrain  roughness  as  the  sample 
spacing  increases,  as  shown  in  Figure  k*     A  filtering  action  takes 
place  so  that  smoothing  of  the  terrain  occurs;  statistics,  such  as  K, 
derived  from  the  triangle  orientations  reflect  this  smoothing. 

Photographic  enlargements  up  to  four  diameters  of  several  test 
sites  were  produced,  sampled,  and  analyzed  using  program  VECTOR.   Since 
all  these  data  sets  were  rectangular  arrays,  it  was  also  possible  to 
analyze  alternate  columns  and  rows.   These  procedures  allowed  the 
examination  of  the  effects  of  grid  size  (sample  spacing)  on  results 
of  program  VECTOR  analysis. 

Testing  of  these  comparisons  for  significance  requires  the  use  of 
non-parametric  statistics.   Strahler  ( 12,14)  showed  that  slopes  have 
normal  distributions.  In  contrast,  the  distributions  of  dips  of  the 
triangular  planes  are  not  normal,  but  vary  widely  from  skewed  to 
almost  exponential  distributions.   Similar  uncertainties  surround  the 
distributions  of  K  and  of  strikes  cf  these  planes.   Transformation  of 
these  distributions  to  more  nearly  normal  form  would  allow  the  use  of 
parametric  statistics.   The  Walsh  test,  described  by  Siegel  (10),  is  a 
suitable  non-parametric  test,  and  was  used  to  evaluate  K  values  measured 
for  Indiana  sites.  The  test  showed  that  samples  collected  on  a  500  foot 
spacing  gave  significantly  smoother  K  values  (at  the  95$  confidence  level) 
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than  did  samples  collected  on  a  125  foot  spacing .   The  investigator 
must  consider  the  effect  of  sample  spacing  on  roughness  estimators 
when  determining  the  roughness  of  any  area  and  set  his  grid  size 
accordingly. 

Most  terrain  types  have  a  characteristic  slope  length.   Statistics 
derived  from  grid  spacings  equal  to  or  less  than  the  equivalent  map 
distance  will  probably  show  much  less  dependence  on  grid  size  than 
will  values  obtained  from  larger  grids.  One  method  of  selecting  a 
grid  size  would  be  to  estimate  the  characteristic  slope  lengths  for  the 
terrain  to  be  studied  and  use  a  grid  size  slightly  smaller  than  the 
equivalent  map  distance. 

Another  method  is  possible  if  the  primary  purpose  of  determining 
the  roughness  is  to  compare  areas  so  as  to  design  suitable  stratified 
sampling  procedures,  as  described  in  a  following  section.  In  such  a 
case  a  grid  size  equal  to  the  minimum  sampling  distance  needed  for 
the  purposes  of  analysis,  applied  to  ail  subareas,  will  allow  their 
comparison  on  equal  terms. 

Estimation  of  the  Hypsometric  Integral 

A  second  method  of  comparing  two  cr  more  homogenous  areas  is  to 
compare  their  hypsometric  integrals.  The  hypsometric  integral  was 
devised  by  Strahler  (l3)»   It  is  a  dimensicnless  measure  of  subsurface 
volume  of  a  drainage  basin.   Strahler  suggested  that,  under  certain 
conditions,  this  integral  would  provide  a  quantitative  expression  of 
stage  of  basin  denudation  and  Sehumm  (9)  successfully  used  it  for 
this  purpose. 
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Strahler's  method  cf  measuring  the  hypsometric  integral  is  somewhat 
tedious  and  this  may  explain  why  this  mcrphometric  measure  has  been  little 
used.   Chorley  and  Morley  (2)  suggested  a  simplification.   They  approxi- 
mated drainage  basins  by  a  regular  geometrical  forms,  a  lemniscate,  and 
were  able  to  more  easily  calculate  the  hypsometric  integral.  However, 
their  technique  resulted  in  values  which  varied  from  Strahler's  values. 
They  suggested  a  transformation  to  correct  their  values. 

Description  of  Program  GCON 

The  availability  of  a  computer  program  GCON,  originally  written  in 
the  MAD  language  by  Professor  W.  R.  Tobler,  University  of  Michigan, 
Ann  Arbor,  suggested  a  second  method  of  estimating  the  hypsometric 
integral.   This  program  accepts  gridded  elevation  values  and  performs 
linear  interpolation  on  them  to  produce  contour  maps  on  the  normal 
computer  line  printer.  Various  scaling  and  contour  interval  specifica- 
tion options  are  available. 

The  program  defines  a  map  as  an  array  of  printing  positions ;  N 
characters  per  line  wide  and  M  lines  long.   Thus  there  are  N  x  M 
character  spaces  in  the  entire  map.  Each  character  position  is  evaluated 
in  turn  and  the  appropriate  interval  list  is  incremented  by  one.   Thus 
a  frequency  table  is  readily  produced  of  percent  area  within  each 
contour  interval. 

Accuracy  of  this  integration  procedure  is  controlled  by  the  fixed 
size  of  the  unit  incremental  area9  a  single  character.   Thus  larger 
map  sizes  will  increase  the  accuracy  of  the  area  determinations. 
Visual  comparison  of  the  computer-developed  contour  map  with  the  original 


Ik 

topographic  map  allows  the  investigator  to  estimate  the  reliability 
of  his  estimated  hypsometric  integral. 

While  the  program  requires  rectangular  map  areas,  analysis  of 
irregular  watersheds  can  be  accomplished  by  setting  all  areas  outside 
the  watershed  to  arbitrarily  very  high  values.   The  maximum  and  minimum 
contour  levels  corresponding  to  the  highest  and  lowest  values  for  the 
watershed  are  included  as  data.   The  program  calculates  the  elevation 
range  and  excludes  all  elevations  lying  more  than  one  range  above  the 
watershed's  highest  point,  or  one  range  below  the  lowest  point.   The 
frequency  table  will  accordingly  show  the  number  of  character  spaces 
within  each  contour  level  in  the  basin.   These  values  can  be  quickly 
transformed  for  plotting  the  hypsometric  curve.  A  new  subroutine 
has  been  developed  to  allow  the  computer  to  make  this  transformation, 
plot  the  curve,  and  estimate  the  area  below  the  curve  -  the  hypsometric 
integral. 

Results  of  Analysis 

The  above  concepts  were  tested  by  re-analyzing  seven  drainage 
basins  whose  hypsometric  integrals  were  measured  by  Strahler  (13). 
figure  5A  is  a  copy  of  Strahler' s  map  of  one  of  these  basins;  Figure 
5B  is  a  copy  of  a  computer-generated  contour  map  of  the  same  basin. 
Figure  6  is  a  plot  of  the  hypsometric  curve  for  the  same  basin. 
Strahler 's  values  are  also  plotted  shewing  the  close  agreement  between 
the  two  methods.   Table  1  tabulates  the  results  for  all  seven  drainage 
basins.   Choriey  and  Morley's  uncorrected  values  for  the  same  basins 
arc  also  listed  for  comparison. 
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FIGURE  5- DRAINAGE    BASIN,  VERDUGO    HILLS,  CALIFORNIA. 
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FIGURE  6-  HYPSOMETRIC   CURVE,  DRAINAGE   BASIN    IN 
VERDUGO    HILLS    (see     Figure  5) 
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TABLE  1 
Comparison  of  Hypsometric  Integrals 


17 


Values   from  Strahler   (13) 

t - —   

Computer 

Estimated 

Hypsometric 

Integral 

• 1 

Approximated 
Hypsometric 

Integral 
from  Chorley 
and  Morley  (2) 

Figure 
Number 

Single-Basin 

Hypsometric 

Integral 

Composite  V 

alues 

Mean 
Hypsometric 
Integral 

.....      ,      . 

Estimated 
Population 

Standard 
Deviation 
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79-% 
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METHODS  OF  STRATIFIED  SAMPLING  OF  MAP  SHEETS 


Most  topographic  quadrangle  maps  cover  a  variety  of  terrain  types . 
Such  areas  may  be  termed  "compound"  as  opposed  to  "simple"  terrains 
which  show  only  one  type  of  landform.  In  statistical  terminology  such 
compound  map  areas  contain  more  than  one  population . 

Suppose  a  researcher  wishes  to  digitize  a  map  area  as  shown  in 
Figure  7A.  He  might  partition  the  map  into  three  "simple"  areas  as 
shown  in  Figure  7B.   The  investigator  may  decide  on  a  total  number  of 
sample  locations  for  the  entire  map,  or  he  may  decide  on  some  sampling 
density  for  any  one  of  the  simple  areas .   In  either  case,  if  he  wishes 
to  sample  all  areas  with  equal  reliability;,  he  must  distribute  the 
number  of  sample  points  for  each  subarea  according  to  its  relative  area 
and  its  relative  variability. 

The  percent  of  the  total  map  area  contained  in  any  subarea  is 
easily  estimated.   It  can  be  obtained  by  a  planimeter  if  greater  accuracy 
is  desired.   Thus,  the  problem  resolves  into  determining  a  suitable 
numerical  descriptor  of  the  variability  of  the  various  terrain  types. 

Terrain  Variability  Measures 

A  perfectly  smooth,  planar,  horizontal  surface  can  be  adequately 
described  by  a  single  elevation  value;  any  planar  inclined  surface 
requires  three  elevation  values,  and  associated  X  and  Y  coordinates,  to 
describe  it.   Gently  curved  surfaces  require  still  more  data  points  to 
satisfactorily  describe  their  overall  orientation.   The  overall  orien- 
tation of  subareas  might  be  described  in  a  number  of  ways.   The  range 
of  the  elevations  is  the  simplest  value  to  obtain,  requiring  no  prior 
sampling  procedure. 
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A.  BLOCK  DIAGRAM   OF    MAP   AREA      (after  Lobeek 


AREA 

TABLE 

SUBAREA 

%    AREA 

1.   bedrock  hills 

47 

2.  ground  moraine 

57 

3,  flood  plain 

36 

B.  MAP  SUBDIVIDED  FOR   STRATIFIED   SAMPLING 


FIGURE  7- STRATIFIED  SAMPLING  PROCEDURES 


TABLE  2 
Results  of  Program  VECTOR  Analyses 
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Test  Site 
Number 

Terrain  Type 

______ — 

Elevation 
Range  (ft) 

Roughness 
K 

log  K 

Varia- 
bility 
v  (ft) 

A.   Turner  Test  Sites 

02393-1 

escarpment 

315 

85 

1.929 

163.5 

02393-3 

escarpment 

278 

82 

1.911* 

145-1 

03022-1 

hills 

340 

397 

2o599 

131.0 

02487-1 

hills 

277 

220 

2.342 

118.5 

02039-2 

hills 

212 

158 

2.199 

96.5 

W-l 

drumlins 

183 

222 

2.346 

78.3 

02039-1 

karst  plain 

166 

723 

2.859 

57.2 

02039-3 

karst  plain 

170 

1329 

3.121+ 

54.6 

02393 

plateau 

81 

1759 

3. 21+5 

25.0 

02393-2 

ridge  moraine  and  dunes 

53 

2892 

3.^61 

15.3 

02329-2 

ground  moraine 

60 

24661+ 

4.392 

13.7 

02329-1 

outwash  plain 

28 

18024 

4.256 

6.6 

02625-1 

ground  moraine 

2b 

30044 

4.478 

5.8 

B .  Stone 

and  Dugundji  Sites  (subar 

ea  numbers  a 

re  figure  n 

umbers ) 

M-29 

Boulder-free  wadi 

25.0 

31 

1.491 

16.8 

M-57 

Pleistocene  lake  terrace 

10.8 

51 

1.708 

6.3 

M-33 

Mi  cr 0 -b  adl ands 

9*0 

117 

2.068 

4.4 

M-45 

Sand-sheet 

10.5 

309 

2.490 

4.2 

U-k9 

Complex  dunes 

8.7 

l6l 

2.207 

3.9 

M-53 

Turret  dunes 

7*0 

133 

2.124 

3-3 

M-4l 

Playa  drainage  channels 

5.3 

233 

2.367 

2.2 

M-37 

Floodplain  mounds 

3  =  0 

269 

2.430 

1.2 

M-39 

Floodplain  ridges 

2.3 

189 

2.276 

1.0 

M-35 

Salt  ploygons 

l.k 

102 

2.009 

0.7 

M-VT 

"Devil's  Cornfield" 

1.15 

975 

2.989 

0.4 

1  M-55 

Playa 

0.47 

374773 

5.573 

0.084 
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A  particular  subarea' s  roughness  will  also  affect  the  number  of 
sampling  points  required  to  describe  it.  It  has  already  been  shown 
that  K,  Fisher's  dispersion  factor s  is  an  ideal  numerical  descriptor 
of  roughness.   Since  K  values  can  range  from  one  to  extremely  large 
numbers  the  variability  factor,  v,  is  defined  as 

(Range  of  elevations  (feet;) 

V  =  "*— S : I ; — ~*- 

log  (K) 

v  can  range  from  zero  for  completely  flat  surfaces  which  must  by 

definition  be  smooth,  to  the  range  in  elevations  for  extremely  rough 

areas  (K  -  0,  thus  log  (K)  =  1).   Table  2  lists  values  of  elevation 

range,  K,  log  (K),  and  the  resulting  v  for  all  twenty-five  test  sites. 

If  the  variability  factor  and  the  percent  area  belonging  to  each 

subarea  on  a  map  sheet,  such  as  Figure  7B,  are  known  and  the  total 

number  of  samples  to  be  distributed  on  the  map  is  decided  on8  the 

following  equation  will  distribute  the.  points  among  the  various  subareas 

v.a.K 

" i  ~sv.  a. 
i  l 

where  i  is  the  subarea  number,  v.  is  the  variability  of  the  ith  sub- 
area,  a.  is  the  percent  area  of  the  ith  subarea,  and  M  is  the  total 
number  of  samples  on  the  map.   Such  a  distribution  of  sample  points 
among  the  various  subareas,  will  result  in  approximately  equal 
sampling  reliability  over  an  entire  map  area. 

CONCLUSIONS 

Homogenous  terrain  types  can  be  compared  using  measures  of  roughness, 
variability,  and  the  hypsometric  Integral.   The  digital  computer  can  aid 
In  the  calculation  of  these  measures.  In  addition  to  their  geomorphic 
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applications  these  measures  could  be  used  as  predictors  of  several 
engineering  construction  costs ,  such  as  highway  grading  costs ,  since 
they  are  associated  with  terrain  geometry. 

The  variability  factor  can  be  used  to  determine  appropriate 
irregularly  spaced  sampling  distributions  of  larger  compound  areas 
covering  one  or  many  map  sheets  .  The  importance  of  these  techniques 
should  not  be  underestimated „  Mclntyre  (7)  states  that  irregularly 
spaced  random  samples  may  better  digitize  topographic  maps  than  gridded 
data  under  at  least  some  conditions .   Computer  analysis  of  these 
compound  areas  by  techniques  such  as  trend-surface  analysis  will  open 
the  way  to  a  number  of  interesting  studies  . 

Map  comparison  procedures,  as  suggested  by  Miller  (8),  are  one 
obvious  application o  Geometric  analysis  of  at  least  more  regular 
terrain  types ,  such  as  drumlinized  topography,  by  the  methods  developed 
by  Whitten  (l6)  and  Loudon  (6)  could  have  useful  applications  to  glacial 
geology,  geomorphology,  and  civil  engineering*  The  author  is  currently 
examining  the  applicability  of  trend  surface  analysis  to  the  highway 
location  problem. 
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APPENDIX    A 
Maps  of  Turner's  Test  Sites 
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FIG.  A-l.     Test  site  02393-1;  Escarpment,  Medora  Quadrangle  (I  24.000) 
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FIG.  A-2.    Test  site  02393-3;  Escarpment,  Medora  Quadrangle  (1 :  24.000) 
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FIG.  A-3.    Test  site  03022-1;  Hills,  Oolitic  Quadrangle  (1-62.500) 
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FIG.A-4.   Test  site  02487-1,  Hills,  Owensburg  Quadrangle  (1.24.000) 
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FIG.A-5.    Test  site  02039-2;  Hills,  Bedford  West  Quadrangle  (1:24.000) 


FIG.  A-6.   Test  site  W-l;  Drumlins,  Weedsport  Quadrangle  NY   (I   24.000) 
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FIG.  A -7.    Test  site  02039-1-  Karst  Plain,  Bedford  West  Quadrangle  ( 1 :  24.000) 
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FIG.  A-8     Test  site  02039-3;  Karst  Plain,  Bedford  West  Quadrangle  (1:24  000) 
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FIG.  A-9.  Test  site  02393-4;  Plateau,  Medora  Quadrangle  (1.24.000) 


Fig.  A- 10.  Test  site  02393-2;  Ridge  Moraine  and  Dunes,  Medora  Quadrangle  (1:24.000] 


FIG.  A -II.  Test  site  02329-2;  Ground  Moraine,  Lafayette  West  Quadrangle 

(1.24000) 
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FIG.  A-12.  Test  site  02625-1;  Ground  Moroine,  Tipton  Quodrangle  (1.24.000) 


FIG.  A -13.  Test  site  02329-1;  Outwash  Plain,  Lofoyette  West  Quadrangle 

(1:24.000) 


APPENDIX    B 
Computer    Program    VECTOR 


$ID 


ELEVATION     MATRIX 


PROBLEM     CARD 

PROB     (COLS     7  -  80) 


FORMAT    CARD 

FMT     (COLS      7  -  80) 


OUTPUT     OPTION     CARD 

I  PLOT    ,        I   SKIP    ,         I   H  1ST     ,      I    TYPE 


SIZE/  SCALE      CARD 

M  ,      N    ,       XSTEP    ,     YSTEP  ,     DIRN 


TITLE      CARD 

TITL     (COLS       7-  80) 


$DATA 


PROGRAM  VECTOR 
SOURCE  DECK 


FIG.      B-l  PROGRAM      VECTOR 


TYPICAL    DECK    SETUP 


FIG.  B-2.  Listing   of  Program   VECTOR  32 

StD    3402*2*1 0**TURNER*     PROGRAM   VECTOR         ~  ~~ 

SEXECUTE         IBJOR 
BIB JOB 

SlPFTC  MAIM 

CX#****x********  ****** W#W#*&#**W*# *#*"### **¥**« -••- -«.j<.j(. .»..-;..«..-.  jiji  j.  «  "-j'.v.K^jf  44Jtjix^_-<.Jt.^jnm^4m.1nt 

C 

C   PROGRAM  VFCTOR  —  PUROUF  UNIVERSITY  VERSION 

C_   THIS  PROGRAM  JNALVSFq  ^llRCAr^  nmir,HMF<^^ 

C  TWO  TECHNIQUES' ARE  USED — 1)  ELEVATIONS  ARF  ANALYSED  FOR  MEAN  QftMSr 

~~C  VARIANCE  »  STANDARD  DFV'ffiT  I  ON  <  FTC.  2)  GROUPS  OE  THP^r  PO]MT«  ARE  USFD 

C  TO  DEFINE  PLANES  WHOSE  ORIFNTATION  IS  DETERMINF^,  McrA5upc-^  0E  VICTOR 

C  STRENGTH  AND  DISPERSION  APf  u?c^  TO  ESTIMATE  ROUGHNESS  •  Tyoc  j 

C  ANDTVPE  2  TRIANGLFS  CAN  °e  FITTED.  A  VAPI"ILITY  cACTOPi  V  i   TS 

~C  CALCULATED  FROM  THE  RANGE  OF  FLrV4T  I ONP/THF  RjOUGHNFSS  MEASURE*  K  • 
C 

£*******-#********#*******#****##****•>(.**•*■■**■  *##-•(..*.!{..* -a****.**.**,  x****,,.  *****:*.***#*..»**** 
C 


C       THIS  PROGRAM  USrS  THF  ^OLOWING  SUFPOUT  INFC _ 

C  TILTS WHICH  CALCULATES  STRIKES  AMD  DIPS  FRO"'  VICTOR  COSI^'FS. 

C  STATS  —  ",'HICH  CALCULATES  STATISTICAL  M^A^wPFSi 

C WULFF  —  WHICH  GENER^T-S  WULFF  NET  ST^R^OGR A"S  ?r  ^  »rrrn<;. 

~~C  HISTO WHICH  GENERATES  HISTOGRAMS  OF  STR I KF  AND  TIF  0  I ctd I nUT j 0mc , 

C 

C*********************  *******************♦**#  *****jf********^#*^*##*jinnnnnf***»*** 
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C  CARD  CPd^p  IS  

C  A)   CARDS  WHICH  APPEAR  ONC17. 

C     1 )  SPATA 

C    2)  TITLE  CARD  ( T I TL )      COLS  7-80 

~C  A  GENERAL  TITLE  WHICH  WILL  APPEAR  AS  A  HEADING  AT  THE  ',-C!|>'tM( 

C  OF  EACH  DATA  SET. 

C 

C  B)   CARDS  WHICH  ARE  REOUIRFP  AS  PART  OF  FACH  DATA  SFT. ■ 

~C     T)  SIZF/SCALF  CARD  (M»N(  XST^P  ',  YSTEP IDIRN) 

C  M-< COLS  1-3) I IS  THF  mumpfR  OF  ROWSt 

C          N-<C0LS4-6)  I  15  TH^  NUMBER  OF  COLUMNS. 
_C XSTEP-(COLS  11-?0).TS  THF  s  1 7F  OF  COLUMN  spACTmg  tm  fi^watiom  ijnTTS* 

C  YSTEP-<COLS  21-30)»IS  THE  SIZE  OF  ROW  SPACING  IN  ELFVATIOM  UNITS. 

C  01RN-JC0LS  31-40)*  IS  AZIMUTH  IN  D^GR^FS  OF  Y  AXI"=. 

C    2)  OUTPUT  OPTION  CARD   ( I  PLOT  •  I SK I P«  I H I  ST ( 1  J ,  I H I  ST ( 2  )  ♦  I TYPf ) 

C IPLOT-(COL  1HIF  IPt  OT=1  «STR 

C  I  PLOT-  (  COL  1  )  »  IF  IPLOT=l  .STEREOGRAM  OF  MQPMAl  WEC-rooe  IS  pDrmyrpPi, 

C  IS<IP-(COL  5)  I  IE  ISKIP=1  »A  LISTING  OF  TRIANGL^  CTF!^r''.  A»'^  njoc 

C                                        IS  cUPPRFcSed. 
_C IHlST-(COLS  10-1  IMF  IHIST(1)=1.A  HISTOGRAM  OF  Pipe  TS  PLOTTED. 

C    •  IF  IHIST(25=1«A  HISTOGRAM  OF  STRIKES  IS  PLOTTFD. 

C  I  TYPE- (COL  12)  «   IF  ITYPE=1.  TYPE  1  TRIANGLES  ARE  FITTED. 

C                              IF  ITYPE=2t  TYPF  ?  TRIANGLES  ARf  fitted, 
_C IF  ITVPEg3«TVPFS  1  AND  2  TRIANGi  FS  ARF  FTTTFD. 

C         (   IF  THIS  CARD  IS  FLANK  STANDARD  OUTPUT  IS  PRODUCED. TYRE  \     TRI- 

C         ANGLES.A  LISTING  OF  ALL  TRIANGL^  STRIKFS  A"'D  DIPSiAMD  MO 

C         STEREOGRAM  OR  HISTOGRAMS,   ) 

C     3)FORMAT  rAPD  (FMT)    COI  S  7-RO 

~C"  SPECIFIES  FORMAT  OF  DATA  SET  WHICH  FOLLOWS 

C    4)  PROBLEM  CARD  (PROR)      COLS  7-80 

C  A  TITLE  WHICH  WILL  SERVE  TO  IDENTIFY  THE  DATA  SET. 

C    5)  SERIES  OF  D&TA  CARDS  GIVING  ELFVATION  MATRIX* 

C******X****XX*  XXX****  ***********  ******  XX*  XX*  XXXXXX****-"   f-f  ••  ^«-*#-^^-v.#^#^.4f^#1(.^4(.^jj.^Mt 

C 

C   PROGRAM  MODIFIED  FROM  R.D.  HQPSON.  KANSAS  COMPUTER  CONTr  I  3L'T  T  ON  14 

£*************#************#*  ****#****-*##*******-S- •*■*•**•*■  XXXXXf-*?****************** 

_£ A.  KEITH  TURNER ; „ 

<-  PURDUF  UNIVERSITY 

C  AUGUST  1967 

C 

e-^«»«**W*^■^!^-^?-^<■?^^»-?«<■*»K»«»^}»»J-.;«^»:^)p^;»«»y,^^^^gy^y^^■^-»^^>.  «;.;p  *  •+*  jy;  "  1*4  i$xx***#-# 

DIMENSION  E<  20.20)  .XNORMISOO)  .YNORM(ROO)  »2NORM(800J  i<">PSCr")0>  » 
lSTR<a00)»TlTL(l2>«PMT(l^).PROB(12) 

DIMENSION  IHIST<2)      _^_^_ _____ 

<_      -*##*#*  ****-■■  xxxxxxx*  *********  x***#*#**#** 

C   READ  TITLEtMASTERlANO  DATA  FORMAT  CARDS 

C     #**x*xx**x ******* *xx****xx**x* ♦♦##**#*#» 
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READ(5.108>  TITL 

108  F0RMAT<6X« 1  2A6  > 

500    READ  (5.109)  M iN .XSTEP. YSTEP. D IRN 

109  FORMAT  (?I3.4Xt3F10,  ) 

RFAD(S«  1000)  I  PLOT*  !  SK  I  P  •  (  I  H  I  ST  (  I  )  .  I  =  1  »  2  )  t  ITYP!^ 

1000   FCRMATf  I  1  i3Xi  I  1  •  4X t  2  I  1 .  II  ) 

READ  (5.108)  FMT 
C     *****#*********#**#*■*##*#*#*#**#■»###***# 
~C WRITE  OUT  HEADINGS 

C       ********#*********#*#******#***»******## 

WRITE  (6.110) 

110  FQRMATt 1H1 j 10X.30HCOMPUTATIQNS  BY  PROGRAM  VECTOP  ) 

WRITE  (6,H6> 

116  FORMAT  (1  HO  »49H— —  //) 

WRITE <6« 106)  TITL 

WR I TE ( 6 • 1  16) 

WRITE(6» I  OH)  M.N.XSTEP.YSTeP.DIRN.FMT 


100  FORMATdHOt  14HVALUES  READ  IN//.12H  MO  OF  RO<"<5=  «  I  3  .  1  2H  S'O  OF  COLS=« 
1I3.7H  XSTEP=«F8.2.7H  YSTEP=  »  FS.2  ♦  1  8H  Y-AXIS  D  I  P.ECT  I  ON  =  •  F8.  2.  //\  3H 
?DATA  FORMATS »12A6  ) 

WFAD  (Ft  105)  PROB 

'•'RITE<6t 116) 
WRITE C6» 130)  PROP 

130  FQRMAT( 1HC.32HCOMPUTATIQN  R'SULTS  FQR  SugAREa  il?A6  ) 

WRITE  (6.116)  '     

C       ****#**#*************#****#***#*#■*■•**■*■*..!!..» 

C   READ  AND  WRITE  ELEVATION  MATRIX 

r     #**************#*********.*****#*•><■*.><..)(..*#.*.* 

DO  31  1  =  1  «M 

31    RFAD(5|F»'T)   (E(  I  t  J)  .  J=l  .N  ) 
'."RITE  (6*  131  ) 

131  FORMAT( 1H0. 23HQPIGINAL  ELEVATION  DATA  // ) 

DO  32  1=1 «M 

WRITEC6.132)  <E< I iJ) «J»1 »N> 

132  FORMAT<12F10,0> 

32  WRITE(6t 133) 

133  FORMAT(lHO) 

CALL  STATS  (  M«Ni 20 . 20  tE«  EBAR I  EVAR  »ES  I G  t  EMAXi  ="M  I  \< ) 
WR I TE  (6.116) 

WRITE (6.210) 

210  FORMAT(lH0.3BHSTATISTICAL  MEASURES  OF  ELEVATION  DATA  > 
WRITE  <6«133> 

WRITE (6, 208)  EBAR »EVAR. ES IG «FMAX« EMI N 
WRITE  ( 6  »  1  1 6  ) 

"~C        ***************  if*  ♦*#****-)(.->!.*  *****♦##**  #*.)(. 

C   GENERATE  NORMALS  TO  TRIANGLES 

C       ***********-"*************  **#*■*•»■*#•!!■*  **■*■*■* 

'.'MaW-l 

NN»N-1 
C     ****#**************#***********♦**#*##** 
C    I TYPE  DETERMINES  TYPE  OF  TRIANGLE  BEING  FITTED 

C     ******#*-**********■•**■**♦#*■*#****#**♦**■##*     

IF  ( ITYPE.LT.l .OR. ITVPE.GT.3)  tTYPF=l   ' 
GO  TO  (51 .E2.51 ) » ITYPF 
51      K  =  l 

C ******  ***»******»**»*»»»**»»».i»  ■»»»*»*»■»-»* 

"C    DO  LOOP  THDOUGH  50  FITS  TYPE  1  TRIANGLES 
C    (  FIRST  TRIANGLE  IS  DEFINED  BY  POSITIONS  ( 1 * 1 )  .  < 1 ♦ 2 ) ♦ ( 2 . 1  >    ) 

C       **********#***#***#*■*■*#********#**■**■**** 

NNK*1    

''/RITE  (6.555)  NNK 
555    FORMAT  (1H1.32HTHE  PROGRAM  IS  NOW  FITTING  TYPE  lIllllH  TPIlNGLr*») 
WRITE  (6.116) 

DO  50  1=1. MM  

M»I  +  1 


DO  50  J=l »NN 

Jjaj+l 

A=(E( 1«JJ)-E( I j J) )*YSTEP 

B"5T_TT7JT=ETTTTT5  )*xst^p 

C=XSTEP*YSTEP 

r*  I V=SOPT  (  A**?+R**?+c**'  ) 

XN0RM(<)=A/DI V 

YNORM(K)=B/DIV 

_NORM<K)=C/DIV 

K=K+i 

A  =  (E(  I I.JJ )  -E  (  I  I «  J  )  > *YSTEP 

B=(E( I.JJ)-E( II.JJ) )*XSTEP 

C=XSTEP*YSTEP 

D I V=SQRT ( A**2+8**2+C**2 ) 
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XNOPM(K)=A/DIV 

YNORM(K)=B/DIV 

ZNOR^<K)=C/DIV 

50   K=K+1 


NK=K-1 
GO  TO  57 
5?      K=l 

_£ ««<-innt»-:,";-;-t-::-:,-;->-aiHi-;inni-lni-nnn 


C    DO  LOOP  THROUGH  53  FITS  TYPE  2  TRIANGLE? 

C     (  FIRST  TRIANGLE  IS  DEFINED  RY  POSITIONS  <l.l)»(2.1).<r>.?) 

C       ********************  **»*******#**•♦♦*#*»# 

NNK  =  g 

WRTTE(~6T555l     NNK 
WRITE<6« 116) 
"0    53     1  =  1  I  MM 
11=1+1 


DO  53  J=l INN 

JJ=J+1 

A= < F !  I  I . J ) -E (  I  I • J J )  ) #YSTFP 

B=(E!I.J)-E!II.J) )»XSTEP 

C=XSTEP*VSTEP    "!—   ~ 

D  IV  =  SQRT( A**2+B**2+C**2 ) 

XNORM(K)  =A/DIV 

YNOR'.'fO  =B/DIV 


ZNORM(K)=C/DIV 

K=K+1 

A=<F(I,J)-E<I.JJ)  )*YSTFP 

B=!E(  I  I . JJ)-E(  I  •  JJ)  )*XSTEP 


C=X"STEP*VSTEP 

D  IV=SQRT  (  A  ■»■■*?+" **p+c**p  ) 

XNORV(K) =A/DI V 

YNORM(K) =F/DI V 

ZNORM<K)=C/DIV 
K=K+1 
NK=K-1 
57      IF(ISKlP-l)   58.59.58 


r  ***********  ***********************«■*■;:- #-!!--)f 

C   CALCULATE  STRIKF  AND  DIP  Or  ALL  TPIANGL^S 
C   BUT  IF  ISKIP=1  DO  NOT  LIST  VALUFS 


58     WRITE (6.200) 

200    FORMAT  < 1H0.31HSTRIKE  AND  DIP  OF  ALL  TRIANGLFS  // ) 
WRITE (6. 201 ) 
20  1  FORWATllHOi  IX . 35HTR I  ANGLE  NUMBER STRIKF DIP 


CONTINUE 

DO  60  K=l .NK 

AL<=HA  =  XNODM(K) 

EETA=YNCRM!K> 


GAMMA=ZNORM(K) 
CALL  TILTS! ALPHA  . BETA  t GAMMA «  STR I KE  »  D I P  »  D I RN  > 
DPS(fO=DIP 
.  __   STR!K)=STRIKE 

5T3     contTmue 

IF(  ISKIP-1  163.69.68 

68  '"RITF(6»202)  (Ki«TP(K)  |DPS<K)  »<=1  .  NK  ) 

2"g FORMAT; 1H  .3(7X.  1.1.  1  1X.F8.4.2X.FB.4)  ) 

WRITE  (6.115) 

69  CONTINUE 

C       **************************************** 

C   COMPUTE  MFAN  AND  STANDARD  DEVIATIONS  OF  STRlKFS  AND  DIPS 

C      ************************-;:-*************** 

CALL  STATS ( 1 .NK. 1 . 800 . STR • STRBAR . STRVAR . STPS I G I STR MAX . STRM I N ) 
CALL  STATS! 1 » NK . 1 iSOOiDPS . DPSBAR iDPSVAR »OPSS IG.DPSMAX.DPSM I N J 
WRITE (6. 207) 


207  FORMAT! 1H0.44HSTAT1STICAL  MEASURES  OF  STRIKES  OF  TRIANGLES  > 
WR I TE ( 6 • 2  08 ) STRBAR . STRVAR . STRS I G I 3TRMAX • STRM I N 

208  FORMATOH  ,5HMEAN=  »F10.  4  «  1  2H    VAR I ANCE= .F 1 0 .4 » 22H    STANDARD  DEVI 
1  ATI 0N=.F1 0.4.1 1H    MAXIMUM=.F10.4.1  1H MINI  Ml  IMa  .PI  ft  .4  ) 

WRITE (6. 2 09) 

209  FORMAT! 1  HO . 4 1 HSTA TI ST  I C AL  M^ASUR^S  OF  DIPS  OF  TRIANGLFS   ) 
WRITE! 6. POP  > DPSBAR .DPSVAR.DPSSIG.DPSMAX.DPSM IN 

WRITE  (6.116) 

— t=t — 


c        **********************###*##*#*#******** 

C    CALL  HISTO  FOR  HISTOGRAMS  OF  STRIKES  AND  DIPS 

_C *****»*»******»*»*»******♦»**»*»»**»«**■* 

IF!  IHIST!  M.NE.i  )  GO  TO  220 

CALL  HISTO (NK.DPSMAX.O.O. DPS »PROT» I »NNK) 
220     1=1+1 
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IF< IHISTt I ) ,NE. 1 )  GO  TO  221 

CALL  HISTOCNK,STRMAX,STRMIN»  STR.PPORl  I  ,  NM'<) 
221     CONTINUE 
C       *****************************  »***S*t«»Ht 


COMPUTE  "FCTO?  PTPE\'GTH 
TOTX=0.0 

toty=o,o 


TOTZ«0.0 

DO  70  K=l iNK 

TOTX=TOTX+XNOPM  <  K ) 

TOTY=TOTY+YNQRM  ( l<  ) 

TO   T0Tz"  =  T0Tz  +  ZNOR'1(K) 

2R=T0TX*T0TX+T0TY*T0TY+T0TZ*T0TZ 
C     ###*#***«********#******** ************** 

C   VICTOR  STRENGTH  SHOWN  ^y  c  ;  TH~R  Rl  OR  P  *  ":  I  SPrRS  I  ON  r^Y  ~~.CT'< 

""£      ***********-■;-*****************-*********** 
R1=S0RT(ZR> 
FNK««NK 

R=R1/FNK 

E"5TK=(FNK-1.0  5/'<*NK-Dl  ) 

WRITE  (6i=c6) 
Kc,:,    rQDMiT  (  1  HI  ) 

WRIT1^  (6,116) 

WRITE (6,204)  RltR 

204    FORMAT) 1  HO . 35HMF ASURES  OF  VECTOR  STRENGTH  AR^  Rl «»F12 »4 »7H|      R»» 
1F9.6,34H,  WHFRF  R=P1/NUMBER  OF  TRIANC-LFS.     > 

WRITF(6«20g)     FSTK 

2^    FORVATI  1H      i^RHVFCTOP    D  I  SP^R^  I  ON  i  ME  ASURFH     BY     FI^HFR?     nf,prrj5if)M     FAT 
1TOR,  IS»F1 2.4,40H.        LOW    VALUES     INDICATE    VARIA^LF     TERRAIN,  ) 

V=(EMAX-EMIN)/AL0G1  O(ESTK) 

VPITE_(5«300)  V 

~7C~0    FORMAT)  1H0«60HTHF  VARIABILITY  FACTOR,  V  =  FLFVATtCM  RANG<=VLOG  ( ESTK 
1)   ,   IS  »F10»3> 
C       ***********■"**************************** 
C  COMPUTE  VECTOR  MEAN)  AND  MEAN  PLANE  ORIENTATION 
XT — DIP-CTION  COSINES  OF  MEAN  VECTOR 

C       ******#****#**#***#*********-:(■*********** 

ALPHA=T0TX/R1 

R^TA=T0TY/R1  

GAMMA  =  T0tZ/P1  ' 

WRITE< 6*999)  ALPHA | BETA , GAMMA 
999  FORMAT) JH0*?1HTHE  DIRECTION  COSINES  OF  THE  MFAN  VFCTO~  AF1^  ALPHAS! 
lg10»5,  7H   BETA=,F10«5«  8H   G  AMMAs  ,  F  1  0  t5  ) 

CALL     TILTS     (ALRHA,"FTA,".AM"A|.STRll<EiDlP«DlRNl 

WRITE (6*206)     STRIKE* DIP 
205     FORMAT) 1^0 ,5°HTH-     ORIENTATION    Or     THE     MEAN    PL'NF      IS     AS     FOLLOWS--- 
1STRII<E=jF8«4j6Hj     DIP=*F"»4) 

WR I TE  (6*116)  ' 

C       ********  ***********  ********************* 

C    "'HEN  IPLOT=l,  CALL  WULFF  TO  PLOT  A  STEREOGRAM  OF  THF  TRIANGLES 

C     **************************************** 


IF(  IPLOT»NE«l  )-  GO  TO  50l 
CALL  "'ULFF  (  NX  ,  XNORM  •  YNORM*  ZNORM  I  PROP  I  NNK  ) 
501     IF  (  ITYPcr.^Oi3»  AN1.NNK«-Q.  1  )    GO  TO  52 

GO  TO  500 

TM3  ' 

ElRFTC  TILT 

SUBROUTINE  TILTS (ALPHA, BETA .GAMMA , STR 1 KE » D I P ,D I RN ) 

C  *  *********  *  *-***■***■»***■**********»*****♦******»»****»*********»»***»****»»»***»» 

~C       PROGRAM  TO  DETERMINE  STRIKE  ANO  DIP  OF  PLANES 
C  ***************************************-»-»  ****4*  ******************************** 
ANGG=ARCOS(GAMMA) 

AC-SORT (  ( ALPHA ) *»2+ ( PETA ) **2 )  *S  1  GN  (  1  .  0  ,  At  PHA ) *S I GN (1  . Q . BFT A  ) 

COSTH=BETA/AC 

STRIKE=( ( 1.5708-ARCOSJCOSTH) ) #57.2957) +D1 RN 

DIPs  ANGG#57.2957 

RETURN ' 

END 
SISFTC  STAT 

SUBROUTINE  STATS (ROWS* COLS,  I  I  » J J • X i XRAR »XVAR , XS I G «  <f  AX»XMtN) 
C************** ******************************** **»**»*»*»*  ***»»*»****»»»»*»*»»» 

~C  PPOGRAM  TO  FIND  STATISTICAL  MEASURFS  ~ ~~ ~_ " 

C ************************************************** ********  »*«#*«**####**«*#*«#« 
INTEGER  ROWS,  COLS 

DIMENSION  X(  I  1  i  JJ) 

^NOsROWS^COLP  """" 

XSOP=0,0 

XSUMsOiO 
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XMAX=X(1  .  1  ) 
XMIN=X(  1  »  1  ) 
DO  10  1=1 .ROWS 
OQ  10  J=l .COLS 


TF  (X(  I  4  J  )  .GT.xMAX  )  XMAX=X(I»J> 
IF  (X( I t J) .LT.XMIN)  XMIN=X(I«J) 

S5flR=XSQB±XJ!  I  «.H  £X  (  I  ..n 


XRAR=XSUM/FNO 

XVARs>(XSOR-?.0*XRAQ*XSUN'  +  FNO*XQAR*XPAP)/F\JO 

XStG=SGRT(XVAR) 

RFTURN 


THIS  SUBROUTINE  PLOTS  A  SCATTER  DIAGRAM  Oc  A  GROUP  OF  DIRECTIO!" 
5SINES  •   THE  DIAGRAMS  MAY  BE  USED  AS  OVERLAYS  ON  A  HO  CENTlMFTEf 
[AMETER  WULFF  NET  •   THE  NUMBER  IN  EACH  RECTANGLE  I  NO  I  CAT FS  TH~ 
JMP.ER  _OF  POINTS  FALL  I MG  I N  THAT  PFCTANGI.F  . 


ENH 
SIPFTC  '."ULF 

SUPROUTINE  WULFF  (  KOUNT?  «  TX  •  TY  .  TZ  <  PROT  <  I  TYPE  ) 
C****»**********  ^»*****-^****-s<  j:  y- **-■;.-■<.*#.::-•;-.» .;;  f:  ■;:-  *  -;:-  j,- **  ** * * *  „  „ 

C  THIS  SUBROUTINE  PLOTS  A  SCATTER  DIAGRAM  Oc  A  GROUP  OF  DIRECTION 

C       COSINES 
C      Dl 
C_      NUM 
C 
C  ****************************************  *************.):  •  >-:■  •****-:;.*  •.-.■•......■..•.  g..^.**^ .»..». #_ * 

C  ***************  *  *  *  *  *  *  *  *  -x-  *  *  *  -:<■  -v-  *  *  *  -■  • * *  *  *  *  -v.  *  *  * *  *  *  *  ■<.:-  *  *  *  *  -it  **-'•-  *  •'•   --  ••  ^  -.  *  *  -••-  *  *  *  -y  *  *  *  *  *  *  *  *  *  # 

C THIS  SUBROUTINE  M£J2JLElFP  A-T^R  SUSSnUXlNE  '."DLFF  in, 

C       LOUDON. T.V.-1964-PCOMPUTER  ANALYSIS  OF  ORIENTATION  DATA  IN  STRUCTURAL 
C       GEOLOGY<s-OFFICE  OF  NAVAL  RESEARCH » G-OG.  RRANCH.TFCh.   ^FT.13. 
C      ONR  TASK  NUMBER  389-1 35 •   129  PAGES. 

C*»***»*****»  *  *  **************  *********»»*********»»*»*■;<.»»»  ******  *****  *********  *  * 
~C~  A.  KEITH  TURNER! 

C  PURDUE  UNIVERSITY*  1967  • 

C**»»»**********r **************  ********  *  "-ft************  ****>..y.  ;;.;:■>■».  (.;    .,-  »  .-  Ijjj  .^  ^  »» » ..^  jt 

DIMENSION  I  ROW (48 .40) »TX ( 30  0 ) 4 TY ( 800 ) » ft ( 800 ) 

DIMENSION  PROP (125 
C     **************************************** 

C    INITIALIZE*  PUTTING  -9  IN  THE  AREA  OUTSIDE  THE  CTFpc-QGPAM  . 

C       ******  -'-  *  -"-  *  *  -"-  *  ->-  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  4J  *  *  *  *  *  *  *  *  * 

WR I TE ( 6 • 1 05 ) 
C 
£  *■¥■■%■  -"-  ■*■  -"  # » -%-  -£-  ■#■  -f- -y-  -y-  -fr  ■y- -*■  -,-  -f  -::-  -v-  ';-  ^  -y-  -v-  -k-  -^ -;:-  x-  ->  y-  -jf  ■# -•;■  jc  -y-  -^  -y.  -x-  -y- 

~C   SEE  COMMENT  ON  CONSTANTS  BELOW 
C     **************************************** 


c 


P=l. 0/(3. 937*6.0) 


0=1 .0/(3.037*5.0 ) 
C   ************************************************** 

C   FOR  POINTS  ON  THE  C  I  OJMFFRrNCc  OF  TH^  PR  I  M  I  T  I  V="  CIRCLE  I 

_C (MP/?)**2+(NQ/?)**?  =  1 t     BY  THE  THOREM  OF  PYTHAGORAS  . 

C    **************************************************** 
DO  6   M«l .48 
R=M 

R  =  ARS(R-?4.0) , 

DO  2  N=l .40 
2    IROW(MiN)*0 

N?  =  SORT < ( 1 .05-(R*P>**?> /0**2)  +  0.5 

N3  =  ?0  -MP 

N2  ■  20+N2 
DO  7   N*l »M3 
7    IR0W<M»N)a-9 

3 DO  6^N=N2»40 

S IR0W(M»N)  =  -9 
C    ************************************************ 
C      COMPUTE  THE  NUMRER  OF  POINTS  IN  EACH  RECTANGLE  OF  TH~  STEREOGRAM! 

C ******************************************************      :-*■(•** 

DT5~3         I  =  i«k:0UKjT2 

UP  »TZ(t) 

IF(UP)  1.28.28         s 

1    TX<I)»-TX<t) 

TVI!)»-Ty(I) 


c  ************************************************** *******#**^**^f 

C      IN  THE  NEXT  TWO  EXPRESSIONS.  3.937  IS  THE  RADIUS  OF  THE  NET  IN  IM^ES. 

£      6*0  IS  THE  NU^FJEP  OF  LINE  OF  OUTPUT  PER  INCH.  5.0  IS  TH^  NUMFFP  n' 
~C     TWO-DIGIT  CHARACTERS  PER  INCH*  0.5  IS  A  ROUNDING  FACTOR. 

C  *****************************************#*******#*****##*#**.# 

28       M*<  <TY( 1 >/< 1.0  +  ABS(TZ< I  ) ) ) > *3.937* ( -6. 0 ) )  + I .0+3. 937*6 . 
_^ Na((TX(l  )/<1.0-f-ABS(TZ(I  )  )  )  )  *3. 937*  (-5.0  )  )  ±1  .0+3.  937*5. 

C 

C  ************  ****** ******jf^*>*********##^M(.* 

C  DOUBLE    THE    DISPERSION 
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C  ***#******■**#***#*****■**■**#  #**##*#****  ***-j 

c 

IF(M.GE.2F)     GO     TO    3  1 

M=2*M-24 


IF(N.LF.20)     N=2*N-20 
IF(N.GE.2  1)        N=2*N-21 
GO    TO     32 

IFf'N.LE.20)    N=2*N-?0 

IfJn.GF.21)        n=2*N-21 

32  IFfM.LF. 1  )     M=l 

I F 1 M  .  GF  •  4  8)        M  =  4P 

rFTRTCETI  )       N=l 

IF<N.G^.40)         N  =  40 

3  I  ROW  ( M  •  N  >  =  I  PO1'1  (  M  i  N  )  +  1 

Q  -3f  ■#  3f  "-  -V-  -ft-  ■*£■  -X-  -*£  -K-  -*f  -Jr  -Y-  ■>!■  ■%  -it-  "  'f-  -&  -K  ->  -?r  *  -^  -^-  ■&  ■&  •£  -if  -5J- 

~C        CONVERT     THE     NUMBER    OF     PO I NTS/RF CTAN^LF     TO     PERCENT 

C  ***^#****#***^#***#***#*#**#*##*#*#*#****#*.!j. 

DO    33    M=l (4fl 

00    33    N=l .40  

33  I  ROW? M.N)  =  (  inO.*FLOiT(  lPOW(M,M)  )  /FLOAT  (  KOI  INT?  )  > +0  .  S 
C             ****************************  *******  ** 

C  PRINT    THE    ^CULT? 

_     "r>6  4  W=l«24~  """'"""""' 

4  WRITE (6* 100)   <  I  ROW < M.N)  »N=1  .20)  I  I   I  ROW ( M 4 N > t N  =  2 1 « 40 ) 
WRITE (61 1 04) 

00  5  '-i  =  25  «4P 

=5    WRITE  (6«  100)  (  IPO'"  (MiN)iN=1i?(1]i(  IR0W<WiN)  |N  =  ?1  »40  ) 
'"PITF(6.  1  06) 
WRITE  (6.109)  PROP.ITYPF 

WRITE<6t  101)  KO'JMTp 

WRITE (6. 103) 
WRITE<6i 108) 
WR I TF  (6*107) 
C  

~C  ***************************************-* 

ICO    FORMATMH  t  20X  •  20  I  2  1  1  H,  ,  20  I  2  ) 
101    FOPWATI 1  HO  1  14  1  1 3H  MFASUREMFNTS  ) 
103   fORMAT(lHC«lllHSTFRFOGRAM  SHOEING  DISTRIBUTION  OF  EQlMISj  W    ££HCE 
INT.  PROJECTED  ON  LOWER  HEMISPHERE  OF  20-CFNT I MFTER  WULEF  NET   V 

104  FORMAT  (  19X.P5HWFST. . . . .  (3)  .  ..,.,, 

1 FAST  (2)  ) 

105  FORMAT  < 1H1 .57X.5HN0RTH  ) 

~TO^ FORMAT(55X.HHS0UTH(  1  ) — 5 

107    FORMATt 1H0.38X.45H ) 

108   F0RMAT(  1H0.9RHDISTRIRUT  ION  HAS  REFN  FXAGGFRATFn  ny  DOUBLING  FAST-W 
1EST  AND  NORTH-SOUTH  D  I  STANCES  FROM  THE  CENTER.   ) _______ 

tov       formatc1h0.  hhsubarea   is   «12a6»17h  tr i angle  typf  s   ,11)        '      ~~ 
return 

END 

SIRFTC  HIST 

""5UBROUTTKTE  H 1 5TO ( N « XMAX »XM I N f  X I PROB  » KK *  I TVPF ) 


£***********************#**#  **##*  ##■!!-###.*.##  K^y.  *************         f.^tf  ****  ******#*♦»*##* 

C    THIS  SUBROUTINE  DL0TS  HISTOGRAMS.  fACH  HISTOGRAM  HAS  40  INTERVALS 
C     AND  THE  HIGHEST  CLASS  IS  SET  EQUAL  TO  100  ARBITRARY  UNIT" 

c^***********************************************************^^^^*************** 

DIMENSION  K<40) .X(N> . IHIST(2> »PROB< 12) 

DATA  STAR/lH*/,( IHIST(KK) »KK=1»2)/6H  DIP   »6HSTRIKE/ 

C **************************** ************ 

~C      CALCULATE  SIZE  OF  INTERVALS  «  SFT  NUMBER  OF  OBSERVATIONS  P"d  INTFP''AI   TO  0 
C     **************************************** 
XTERVL=(XMAX+1 .-XMIN)/40. 
DO  20  !=■!  .40 

n?o    «;(n«o 

c    ********##*******■»***##*#*##*##»«*#*#*## 

C    PLACE  EACH  OBSERVATION  IN  APPROPRIATE  INTERVAL 

_C *****************  »»*■)<..*»*.»*.»*  »-»****»*»*** 

DO  60  1=1 «N  '      '   '     "      _ 

M=< (XC I >-XMlN)/XTERVL)+l .0 
6  0     K<M)=K(M)+1 

C **************************************** 

~~C  KIND  HEIGHT  OF  MAXIMUM  INTERVAL  «  KMAX 

C    SET  KMAX  »  100  STARS. AND  CALCULATE  NUMBER  OF  OBSERVATIONS/STAR* 
C    #***♦****************#***##*****#***#*** 

KMAX»K(1  ) 

D0  70  j„j  ,40 — 

70      IF(K( I ) .GT.KMAX)  KMAX«K(I) 
FACK= 100. /FLOAT (KMAX) 
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DO  80 

Mai  1 

140 

80 

K(M)=! 

-ACK*FLOAT(K(M) )+0.5 

FACK= 

1 .O/FACK 

C 

****** 

!f  **************************  ******* 

c 

°f  I  N'T  TITt_F? 

c 

****** 

»•****-"- ** 

*■"■*****  *********■! 

******  ****  * 

WRITE 

C6«  1  1 

1  ~) 

i  h  i  .<=  t  (  kk  : 

1  iPPOT 

110 

FORMAT* 1H: 

I  .  1  3' 

-IHI3TOGRAM  OF  iAfti?flH  V/fll 

1 1^3 

cnp 

9  ip  A  orr  a 

1  1 ?A^  » 

WRITE 

(61 : 

130) 

I  TYPE 

130 

FORMAT  ( IH0« 19HTYPE  OF  TRIANGLE  = 

.11) 

WR I TE 

(6.100)  1 

rACK 

ico 

FORMATdHO.  lQi 

H=-ACH  *  REPRESENTS  .F10.3.14H  0=3F»V  AT  T  Or- 

\Cm//Sj 

C 

****** 

V-  >\-  -"-  -'<-  - 

i-  -/■-.'-  *- 

***************  ******  -;;*  ** * 

c 

PRINT  OUT  HI 

:STOGRAM 

c 

****** 

.«.  Jt  -i'- .». .; 

;-*** 

****•: 

i-**«**->'***,';-*****#-a-#4H 

(•* 

XBEGN 

=  X"lr> 

DO  90 

1=1 #40 

XEND  =  ) 

K^EGN+XTERVL 

- 

NN=K(  I  > 

WRITE 

(6.120)  : 

XBEGN  t XFNO « ( STAR . 

1  J=l  ( 

INN) 

1?0 

FORMA' 

T( 1H9.F1 

?.6i  1 

IX.  1H-.  ix.fi; 

=»6»3X. 

90 

XPEGN- 

=   XEND 

RETUR' 

END 

SDATA 
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FIG.  B-3.  Typical  output  of   program  VECTOR 
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VALUES  HtAU  IN 

10  QF  RyuS-i  2ff  NV>  Uf  CUL5=  £u    XST£P*.   500.01*  YSTEP"   500.00  Y-AAIS  OIHECTION-     0, 

DATA  FOHMATBflOr*. 1 | 


COMPUTAtlUN  *ESULtS  FOR  SUOHREa  *EEUSPUHT  !■   fjRUMLIN  FIELU  NE*  YORK  STATE 


ORIGINAL  ELEVATIuN 


515.       531, 


47«. 

»!». 

»TSi 

4*4. 

4*0, 

4*8. 

441. 

♦*. 

»»T.               4». 

♦63.       478. 


435.       432. 


455.       428, 


510. 


•  50. 


530.  4T1, 

441.  463, 

444.  444, 

418.  411. 

417.  41T. 

42?.  417, 


STATISTICAL    MfcAiuMfcS    Of 


ST4NDAHU    ntvIAUO 


S95.OOO0  MINIMUM*       412.0000 


tNE    PMObWAM     Is     .i. 


IM!AN(iLES. 


STATISTICAL  Mt« 


HU  UtVlATinN" 
*U  UEVlATlON* 


46.0760     HAXIMUMa     90.0001     MINIMUM*   .68.9766 
J. 1316    MAXIMUM'    17.2086    MINIMUM*    -0. 


FIG.  B-3.  cont'd. 
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THE  VARIABILITY  FACTOR.  V  -  ELtVATIuN  HANGE/LOG 1ESTK )  .  IS  .     78.016 

THE  DIRECTION  COSINES  OF  THE  MEAN  VECTOH  ARE  ALPHA-   -0.00052   BETA"    0.0001B   GAMMA- 

THE  ORIENTATION  OF  THE  MEAN  PLANE  IS  AS  FOLLOWS—-   STRIKE*- 18.6209.  DIP"   0.0313 


HlSTOflRAM  OF   OIP 
TYPE  OF  THIANQLF.  ■ 
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722  MEASUREMENTS 
STEREOGRAM  SHOWING  DISTRIBUTION  OF  POINTS.  IN  PERCENT.  PH0JEC7E0  ON  LOWER  HEMISPHERE  OF  20-CENT ImETER  WULFF  NtT 
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APPENDIX    C 
Computer    Progrom    GCON 


41 


,   SPECIAL      CARD 


~\ 


(REQUIRED      ONLY       IF  TOUR  -  I      OR      2 


GRIDDED        ELEVATION 
DATA 


FORMAT    CARD 

FMT      (COLS       I    "    72) 


TITLE     CARD 

TITL       (COLS      1-72' 


MASTER     CONTROL     CARD 

ROWS,     COLS,     CON,    TOUR,    LINES,     LIST,    INCHES    ,     IHYPSO 


$DATA 


PROGRAM      GCON 
SOURCE    DECK 


_  _i 


FIG.      C'l  PROGRAM       GCON    -     TYPICAL     PROGRAM      SETUP 


FIG.  C-2.    Listing  of   Progrom   GCON  .  42 

IIBZZ  2402*2*1  C**TURNBR*   PR  I  NTER/HYPSOMETR  I  C 

5i_Xt  v*  \J  T  t  1  J  «_/U'-> 

SIBJOB  ZFRO 

SIBFTC  MAIN 

r ******************#******^**#*******************4f #*****## #^;  ::  't..-(.jtj<.jf*3ns.j(.4,^  *■*»■**•*■** 

C   PROGRAM  GCON' FORTRAN  IV  VERSION.  (TRANSLATED  FRO"  THP  MAn  I  ".N^UAGE  'T  pijp^i  i^" 

C   UNI  VERS  I  TV  )  THE  PROGRAM  PRODUCTS  CONTOUR  "A"?  of  GQipQcn  "A^A  ON  th-  !  n'c 

~~C   PRINTER.  LINEAR  INTERPOLATION  OCCURS  WITHIN  TH?  GRID  SOUARESi 

C********#****** *****************#**************#******##* **#**♦**#**#*##♦##**## 

C   _  ^     

— C"   THE  MAXr-'U'1  MAP  '•' I  DTH  IS  12.7  INCHESiTHF  '"  I  DTH  0F  TMr"  oADr-p,  ThF  MAX  I  mum  NU"- 
C   PER  Oc  COLUMNS  OE  INPUT  DATA  IS  A  FUNCTION'  OF  TH^  cnrririrr.  map  i.iIOTH* 
C    1^  PS=  (COL-CCLS)  /(CCLS-1  )    (  5cr  DFFINITI0N9  =>e-Low  )   fg  L^sc;  THA"'  ?t     F^JFO 

c      TH- bL ?  "?I"TIVG  .?0? !  TIPMS   '''ould   ncr    "ssi^n^d   T'i  ''a^ih   roiir"'   spa^f    *"d    so   am 

~C   UNSATISFACTORY  MAP  WOX5CS  RESULT .  THUS  THE  PROGRAM  '"ILL  MOT  PLOT  A  MAP  IF  PC  Is"" 
C   LESS  THAN  2  «  AND  ACCORDINGLY  THE  MAXIMUM  NUMBER  OE  COLUMNS  FO°  A  map  OF 
C   MAXIMUM  '-'IOTH  IS  42.  IF  MORE  COLUMNS  ARE  R^OUlR^n  THE  MAP  MUST  Pr  MAD15,  Im 
C SECTIONS  BY  USING  ONLY  PART  OE  THE  QATA  IN  EACH  PASS.  THr  MUMPER  OF  COt  UMNS 

TT SMOULJ  B'-.  THE  SAME  IN  EACH  PVSS  -OR  THC  MSPS  TO  ^F  AT  THF  t;iMr  srALr,  y~  nr, 
C   THIS  IT  may  =>f  NECESSARYi  AND  IS  ADVISABLE.  TO  RFPEAT  SOME  COLUMNS. 

Q  *  »  »  *  *  *  »  *  *  *  a-  »  ---  »  ■*  -■(■  X-  *  "-  »  »  *  *  *  g  *  *  »  .*.  <.  Jt  rr  .-;■  »  ^-  »  »  Jt  3.  ^  _■;.  ^_  g  ,y  .«.  »  3,  .y  3,  .„.  .w.  ^  ^  ^  .„.  ^  ^  ^  .„  .„.  ^  ^  ^  ^  ^  ^  ^  g  ^  #  ^  ^  ^  ^  ^  ^  ^  ^  ^ 

C   CAPO  ORDER  IS  

C   A)  CARDS  WHICH  APPEAR  ONLY  ONCE  — 

C     1 )  SDATA 

"7! 

C   P)  CARDS  WHICH  ARf  nrQUIPcn  FOR  EACH  DATA  Sfj  

C     1)  MASTFR  CONTROL  CARD  (  RO".'S  «  COLS  •  CON  t  TO"R  (  L  I  N|rrS  ,  (_  I  CT  <  I  NCHFS  .  I  HVDf.n  ) 

C THIS  CARD  SERVES  TO  SPECIFY  THF  JYRf  QF  map  AND  ASSOCIATE  OUTPUT  

~~C REQUIRED  FOR  THIS  DATA  SET. ' —      

C  ROWS- (COLS  1-2)  •  IS  THE  NUMBER  OF  ROWS. 

C  COLS- (COLS  3-4) «  IS  THE  NUMBER  OF  COLUMNS. 

C CON-(C0LS  5-6) >  IS  THE  NUMBER  OF  CONTOUR  I MTEFV  M  s  «  "AX  I  MUM  NUMFFP=19. 

~C  !  C0N  =  DATA  PANC-f  DIVID~n  ^Y  r>F^in>cT~  Cr'ITC')^  INT^PVAL  ) 

C  TOUR-(COL  7),  [F  TOU^  =  0«3°OGRAM  USES  MAX  AND  "Im  Fn"  DATA  FOR  u00^^ 

C  AND  LOWER  LIMITC  ^^  CONTOURS. 

C IF  TOUR=  1  >  D^S  IRPD  TO"  AND  "OTTOM  ^^^iT^'l'~g  grfn  |m  fti-.i 

C  SPFCIAL  PROGRAM  CARD  (  <=£•=-  B^LOW)* 

C  IF  TO'JR  =  2«  VARIABLE  CONTOURS  READ  |m  (  (_  O'.f  to  MICH  ), 

C  IN  THIS  CASC  SET  CON  ^QUAL  TO  T"r  NUMBER  OF  VALU^t?  TO 

C pf  READ  IN. 

~Z  IF  TCUR  =  3«C0NT0UR  INTERVAL  B^CO'-'ES  1/2  STANrARD  np'/IATIOM 

C  UNITS.  (  MAX=+3  SIGMAt  MIN=-3  SI^'A  ) 

C          LlN'^S-fCCL  B)tIF  L I NES=0«  CONTOUR  lIn'^c  ADf  PRINTED. 
_C IF  I.  IN"-S=1  «AI  IEBJdAIE (^V^N  VM...W~d) -f"iTni|n  aft>..^c  ri|»Tn, 

C  IF  LINES=2«  COMPLETE  BANDS  AQE  PRIflTED  (ALL  PRINT  °OSMS.) 

C  LIST-(COLP) » IF  LIST*i » INPUT  DATA  IS  LISTED, 

C  IF  LIST=0  OR  BLAN'<«  NO  LIST  PRODUCFD. 

_C INCHES-(COLS  10-13)  >  INCHES  IS  MAP  WIDTH  IN  INCHES.   IF  GREATER  THA'-!  1?,7 

C  OR  nLAM'<f  i  THF  PPor;PA»'  uc^c  \r>, 

C  COL  =  MAP   li'IDJH   !M   CulpTTrn^-lMrur^.,   10.0 

C  IHYPS0-(C0L  lE)tIF  IHYPSO=l »HYPSOMETPIC  AmalYStS  p^n^or^vrn  on  DATA. 

C     2)  TITL^  CARD-(TITL)  COLS  1-7? 

"C^  A  TITLE  WHICH  '''ILL  SERV^  TO  IDENTIFY  THf  DATA  S"T» 

C    3)  FORMAT  CAPD-(FORM)     COLS  1-72 

C       THIS  CARD  SPECIFIFS  THE  FORMAT  OF  THf  GRIDDFn  DATA. 

C    4)  GRIPPED  ELEVATION  DATA 

~~Z  5T  SPECIAL  CAPD(S)     (PFOUIR^D  ONLY  WHFN  TOUR=l  OR  2  ) 

C        If  T0UR=1«  THIS  CARD  SPECIFIES  MINIMUM  AND  mAXImum  "^L^VATIOni  VALUESi 

C        ZMIN  (COLS  1-10)  »ZVAX  (COLS  11-?"). 

C IF  T0UP=2<  THIS  CAP^  SPT.IFI^S  A  FORMAT  (INCOM)  TO  B^  |  icffv  Tm  nr;,n|nf;  ;  V| 

~~£  A  SERIES  OF  IRREGULARLY  SPACED  CONTOUR  LINES.  THF  VALUrS  ARE  ARRANGED 

FROM  LOWEST  TO  HIGHEST  ON  SUBSEQUENT  CARD^,  CON=  NUMBER  OF  LINES* 

q»  »*»»<-*»»»  »*****»*********-*»  **  »***»»»*■»»*  rr***-**»*»*»^ 
~C  PURDUE  UNIVERSITY  VERSION    /      A.  KrlTH  TURNER 

C  JUNE     1°67 

C 

C*********  *  *  *  *  »  *  *  *  *  #  *  *   -- J--  -■  ■  J  •- » -i  y  jf  j£  *  *  »  ,"■  -"- » -"-  » -y-  ■,"■  ---  -y-  >- » --;-  -"-  *  *  jt  -•;-  -",  y.  -y  y.  i:  y.  y  .t  .^  ».  »  j,.  j,.  ^  y_  ^  y.  y  x,.f)HHMlx  -..  ^  ft 

DIMENSION  FORM( 10)  .TITL(12)  »SYMTAB<21 ) .PSYM( 1 3 1)  %  7 < 1 00 , 1 00 ) « 
1  INCON( 10) i DRANG (25) «TABE(?5) »ZPREV( 131 ) • SYM ( 131 ) 

DIMENSION  CUMH(27) 

I  NTEGER  ROWS  »  COLS  t  CON  t  TOUR  |  COLMAX  .  COL  >  psyM  ,oe  ,  GR  I  n  .  T^  ■  IT  .  Tl  Q'»  . 

1  SYM , BLANK I PLUS »  MINUS  «  STAR • SYVTA" 

REAL  MINZ«MAXZi INCHES* I NCON 

D-ATA  (SYMTABt  I  )  «  1  =  1  |21)/1H0|  1H1  ,  1H?|  1H1|  1H4|  IHT,  Jue,,  jm"?,  1Hn,  lhn, 
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T      lHAtlH^* 1HC»1HD, ]HP,1HF| 1 HG » 1 HH « 1HI i 1 H J* 1 HK/« RL ANK/1 H     /»PLU«V1H  +  / 
2     »MINUS/1H-/*STAR/1H*/ 
C  **************************************** 

C  Rf^An    CONTROL    AND    TITLE    CAPOS 

~C *  *  *  *  *  *  *  *"  *  *"  *  *  ■'.'■  *  *  -:;-  *  *  *  *  -::■  *  *  5?  -"-  *  *  ************* 

1       READ (5* 1 000)  ROWS  *  COLS  »  CON  •  TOUR  «  L  I  i\II7"S  *L  I  ST  t  ImCHFS,  ihy"SP 
1000   FORMATfT] 2,3  11 (FiiOtlXi II ) 

Or^AD     (5*1010)     TITI 

1010       FORMAY     (12A61  

C       #*****■**■■****  ******  *  ****************  «#**.»!. 

C    ZERO  COUNTERS  AND  ARRAYS  USED  IN  FREQUENCY  TABL*7 

r      »**  ******  -x-  -x-  *  -st  *  *  *  *  -■:-  *  *  *  *  *  -;:-  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  * 

JJ=TT 
MM  =  0 
SI  )M  =  0  §0 

TLO'"  =  0 

TouT  =  0  

DO  5  1=1 • 25 
5       TAEE(I)=0.0 

N=R0WS*C0LS 

~Z  *"*"*"*"  *  *  *""*"*"  *  *~  * *  *"  *  *"*  *  *  *  *  *  *  *****  *  *  *  *  *  *  *  *  *  *  *  * 

C       N  =  NUMBER  OF  DATA  POINTS 
C    R^An  DATA  pop^AT  CARD  AND  DATA 

_C *  *  *  *  *  »  -'t  *  "-  *  *  -■:  *  *  *  -':-  -■:-  -■■■  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  ;t  *  *  * 

READ  (5*10?0>  FORM  "     ~ 

1020   FORMAT  (10A6) 
DO  7  1=1 .ROWS 

7 READ  (SaFORM)  ( 2 ( I  *  J ) * J  =  l *  COLS ) 

"C     *********■>-■  *"  -5-  -"-  *  4'-  *  ■•:-  ■--  -st  *  *  *  *  *  *  *■•**■**##**#***** 
C     IE  LIST  =1  LIST  INPUT  DATA t OTHERWISE  GO  TO  20 
C     ***********  ****************#*#********** 

IE  (LIST.NE.l)  GO  TO  2  0 

WPITE  (6*1030)  TITL 
1030   FORMAT  (1H1»12A6//) 
WRITE < Si  1025) 

1025   EORMAT(lH0tgqXt22HLISTING  OF  INPUT  DATA.   ) 

DO  10  1=1 t     ROWS 
1C     ''.'RITE  (6*1040)  <7(  I  ,  J)  |  J  =  l  ,COLS> 
1040   FORMAT  ( 1  HO* 1  OF  1 0 »3 • (1 X* 1 0F1 0 »0 )  ) 
20     WRITE  (6*1030)  TITI 


WRITE  (6*1050)  N 
1050   FORMAT  OH  * 1 4 *  1 3H  OBSERVATIONS  // > 

CALL  STATS (ROWS 4 COLS* 100*1 00«Z»ZBAR*XP|ZSIG«MAXZ«MINZ) 
C     **************************************** 


C    DET^P'-IINE  CONTOUR  LINE  ELFVATIONS 

C    TOUP  VALUE  DETERMINES  METHOD  USED  IN  DEFINING  COMTOUDS 

C      *********** ************ ***** ************ 

IF  (TQUR  iGTi  3  iQR.  TOUR  il,„E.  0)  GO  TO  60 

GO  TO  (30*40*50) «TOUR 
30     WRITE  (6*1060) 
1060   FORMAT  (25H  CONTOUR  LEVELS  SPECIFIED) 

RFAD  (5.1070)  7MIN.7MAX 

1C70   FORMAT  (2F10.0) 

GO  TO  70 
4C      PPAD  (5*1020)  INCOM 

' WRITF  (fi.lOP.n) , , . 

1080   FORMAT  <  37H  VARIABLE  CONTOUR  INTERVALS  SPECIFIED) 
READ  (5«INC0N>  (DRANG< t ) *  I »1 ♦ CON ) 
ZMIN=0RANG<1 ) 

ZMAX°DRANG(CON) . , , 

GO  TO  70 
50      ZMIN=Z8AR-3*0*ZSIG 
ZMAX=ZBAR+3*0*ZSIG 

C0N»13 

WRITE  <6*1096> 

1090   FORMAT  <30H  CONTOUR  INTERVAL  ■  Ul/2  SIGMA) 

GO   TO   70 
60  WRITE    (6*1100) 


1 100 FORMAT     (26H    CONTOUR    LEVELS    CALCULATED) 
ZMlN=MINZ 
ZMAX=MAXZ 

70  WRITE     (6*1110)     F0RM*R0WS*C0LS*ZBAR*2SIG*MAXZ*M1NZ 

HIO  ^WMAT  <13H  SKIDDED  DATA/^H  LINEAR  INTERPOLATION  WlfHtN  «Rtn  SCU 
1ARES  /«17H  OBSERVED  VALUES  /«1H  10A6/.6H  ROWS  «  iI3«10H  CClls  *  i 
2i3«/t8H   ZBAR   *N»P10»3#10H         ZS1G    ■    *F10.3»10H         MAXZ    «    «FlO«3«10H 

3   MIN2  ■  »F10«3) 

"C    *****************  *************  *********** 
C    CALCULATE  MAP  WIDTH  IN  CHARACTERS  (COL) 
C     ***************** ****#*******#********** 


IF  MrH^I    r$TA    A2!7'     C0L"127 

gJ5P«-grS^yDtIWCHe>,Ue'«'T>     COL.INCHBWIO.C 
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1120 


FORMAT^ 25H  TOO  MANY  COLUMNS  IN  GRID) 


??hv^af!s*;co*s:'>**"' ******** 


"1  130 


11  40 


GRID=PS+1 

FC0N=C0N 

Cl=RANG?/'FC0Kt " 

If«MST'?Ji'J3iNr:I^fg:x;f..ep.o.xiNc.y.Nc  ,CON,PS.  C0L-« 


2)     GO    TO    100 


IF     (TOUR.FOc 
JJJ=C0N+1 

-jkANG{I)n2MtN+FL0ATn-n*Cf' 
RANG2=a,0*RANGZ  I>*CI 

RL0W=ZMAX-RANG2 
RtHGH  =  ZM  I N+R  ANG? 


I5t 


TTTC- 


1  180 


1190 

1200 


^tf^raSfe'^TTaSV,!^    »«    AS   FOLLOWS)  ~ 

Jg^VgHO^^VEN    SYMBOLS    PRINTED     IN    .OOy'  OF    MAP,      , 

IF    <LINES»NE.0)     GO    TO     110 

WRITE    (6*1200)  ' 

JLOPMAT    I^Xfi^tAXtlOHAgoVE    HIGH) 


I  PhU 


vAlIbI  .0-YtNC    ^ ■ 

GRID-C0LS-1 
C0L«3 

N=PS+1 

DO  140  J»1«GRID 


130 


140 


1240 


00  130  !-  =  2»N 
SYM< COL )» BLANK 
C0L»C0L'f  1 

bVM( COL) « PLUS — 

C0L»C0L+1 
N=C0LMAX+2 
SYMf2)«PI  !)■: 

S^SI^-blank 

XSTOP«ltO-XlNC 

^iT(fi,Hi:?29i?r(n"! 


1  IN) 
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C     *¥*¥*****¥**«W»***5Mf  **#*  ***•***#♦*•******* 
C    ESTIMATE  VALUE  OF  ALL  PRINTING  POSITIONS  ON  MAP 
C     *•**************************»***##****#** 
_C LOOP  THPOUGH  VAP  LENGTH  i  ROW  RY  ROW 


*****•»<■*  *  ***  *  *  *  -"-**  *  *  *  *  *  **  *  ***  **  *  -if  *  - 

O      74"      '-".IsJiOnHC 

=MM-1 

IF     (v-'.G^.^O'"S)     YSTOPsl.n 


T55  YP=0.0 

150         C0L=1 
K=l 
IF     (YP.NE.0.0)     GO    TO     160 


DSY"( 1 ) =PLUS 
PSYM(N) =PLUS 

go  to  l7n 
PSYM(  1  >=^LAM'< 


PSY«(N)="LANK 
170    IF  ( <M.EQ.1>.AND«<YP.EQ«O»0))  GO  To  ISO 

IF  ( (MM. GE. ROWS) .AND. (YP.GT.VAL1 ) )  go  TO  ! BO 
I  NO EX  =  0 


GO  TO  190 
180     INDEX=1 

PSYM( 1 ) =DLUC 
p^Y''  (  N  )  =  PLU.C 


C       *****************  ****************  ******* 

C  LOOP    THROUGH    MAP    WIDTH, COLUMN    BY    COLUMN 

C  **********************************»**#** 

190  IIaCOL^-1 

^C    330    J=l  .  I  I 


JJ=J+1 

IF     (JJ.G^.rOLS)     XSTO°=1.0 

IF     UiGTill     71  =  7? 


IF     (J.LE.l)     Z|aVP*<Z(MMi J)-Z<M« J) >+Z(M» J) 
Z?=YP* ( Z ( MM •  J J ) -Z < M t J J )  ) +7 ( M «  J J ) 

ZDELT=Z2-71 

XP=0.0 


195    COL=COL+l 

r  #*#**##*******#**#***■##**####*#*###■*■**#* 

C       CHFC<  FOR  FDGE  OF  MAP 
J* ***:■■**  :  -"  -■•  ■■■--■■--■■.-■•-*.■,•-■■■  -- ^----r- >•--■•- y- y- "-*-;■-**»*»»******* 


IF  (LINES. EO, 2)  GO  TO  200 

IF  (C0L.NE.2.ANC.C0L.LF.C0LMAX.AND. INHFX.NF. 1)   LIST=0 
200 ZEST=XP*ZDFLT+Z1 


IF  (C0L.GT.2)  <=0 
*****-■•- ******************-».->  J,.**  ****.&.* -a.*  .«..&.£ 

C      ELIMINATE  OUT  OF  RANGE  VALUES 

_£ -■:-  *  *  *  *  *  *  *  *  ■>•-  *  -"-  *  *  *  f  »  *  *  *  *  * J  --«-**********»****  * 


IF  (ZEST  .LT.  RHIGH  .AND.  ZEST  .GT.  PLOW)  GO  TO  210 

TOUT=TOUT+l 

IF     (LIST.EO.O)     PSYM(C0L)=3LANK 

IF     (LIST.NE.O)     PSVMtCOt    1=STAR 


GO    TO     32  0 

C  ******-"********************************* 

C     ASSIGN  ELEVATION  ESTIMATE  TO  APPROPRIATE  CONTOUR  INTERVAL 
_C *  g  -"-  -"■  "-  -"-  "-  '•■  ---  *  *  *  -"-  *  *  -"- *  -"-  *  »  ******************** 


"2H5    SUMsSUM+1 .0 

IF  (ZEST.GE.ZMIN)  GO  TO  220 
I=-l 
TL0W=TL0W+1 


GO  TO  270 
220     IF  (Z^ST.LE.ZMAY)  GO  TO  '30 
I=C0N+1 
T&aEJ  I  )=TAgc-<  T  1  +  1  . 


GO  TO  270 
230    IE  (T0UP.E0.2)  GO  TO  2A0 
I  =  (  < ZEST-ZM I N  > /C I >  +  1  .  0 
GO  TQ  26^ 


240    DO  250  1  =  1  «CON 

IE  (ZEST.LT.^PANG ( I ) )  GO  T0  260 
250    CONTINUE 

260 CONTINUE , , 

265 TABE(  I)=TAPF(  I  )  +  i  , 

C     ***#***********♦*******+**#****#**»***** 

C    ASSIGN  PRINTING  CH4RACTFP  TO  EACH  PRINTING  POSITION 

§******»****  ^-***  ******»»**»»****»»»»*»**» 
78     IF  (LIcT.fo.O)  GO  TO  2^0 

IF  (I.LT.O)  PSYM(COL)*MINUR 
Ir  (I.GT.CON!)  PSYM(COL) =PLUS 
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IF  (  I  .GE.O.OR.  I  .LE.CON)  PSYM  <  COL  J  «SYMTAR  (  II 

GO  TO  320 
C      -I;-*************************************** 
C    PRINT  MAP  ACCORDING  TO  LIMES  SPECIFICATION' 


C      **************************************** 
280     IF  (LlNES.^O.O)  GO  TO  ?<-"? 
C     *****-jh:- ******************************  ■ 
C      PRINT  EVEN  .^Y'^^L5 

J***********************************.)"' 


K=I/2 

IF  <  (K*2)  .NE.  I  •  AND.ZEST.LT.ZMAX)  PSYM  ( <*OL  )  "SYMTAS  (  I  ) 

IF  ( (K*2).FQ.I.0R.ZFST.GF.?MAX)  PSYM ( mil =BJ  AMK 


GO  TO  320 

C      ***************************************** 
C       PRINT  CONTOUR  LINES 

_C ****»***»***********»*****************-••- 

"2^5^    IF  (K.EQ.O)  GO  TO  30 
PSYM ( COL )=SYMTAB( I ) 
GO  TO  320 

C ******************»*»***»***»***»»*»*»»* 

~~C  ",'HICH  CONTOUR  PASSES  THROUGH  PPINTING  APC4 

C       ***********************************  x*-%-)!~* 
300    DO  310  1=1 .CON 

ZSTEP=DRANG(  I  ) 

5W!»1 . 


SW2=  1  • 
SW3»1 • 

SW4°1 I 

"TF  (Z5TFP.GE.ZLAST.ORiZSTEP.GE.ZPRFV(C0L>  )  S>"1=0. 
IF  (ZSTEP.LE.ZEST)  SW2=0. 
IF  (ZSTEP.LT.ZLAST.OP.ZSTEP.LT.ZPP^VfCOL) )  S»'3  =  0. 

} Z   S  Z^.TE- * r'_ ' 7-\ SZ \  _ g.- A B n  * 

IFC (SW1  .EQ.O..AND.SW2.FO.0.  )  .OR.  (  SW3.EQ.0  •  •  AND.  SV»4  ,FC>»0  .  >  )  GOTO  310 


PSYM ( COL )=BLANK 

GO  TO  320 
310    PSYM(COL)«SYMTAB( I) 
"32D ZLAST  =  7EST 

Z"REV(CCL ) =  ZEST 

XP=XP+XlNC 

IF  (XP.LE.XSTOP)  GO  TO  195 

~3"30 CONTINUE ~ ~ " 

WRITE  (6.1240)  (PSYM(L ) .L=l «N) 

XSTOP=XSTOP-XlNC 

JJgQ 

YDeVP+YINC 

IF  (VR.LF.YSTOP)  GO  TO  150 
340    CONTINUE 

.  'r-gHg.  'a'.  Ig^O)   (SYM(L)  »L  =  t  IN) 
"C"     *******************#*S*#*#***4**#*#*#*** 
C    CALCULATE  AND  PRINT  FREQUENCY  TABLE 
C     **********#******#****##**##*####*#***#* 

WRITE  (6  ♦1250)  TOUT 

1H50   FORMAT  ( 16H1FRE0UENCY  T ABLE • /♦ 1 2X . SH I NTERVAL  .16X.5HUNITS  .6X.10HP 
1ERCENTAGE  * 4X • 1 OHCUMULAT I VE  ./.16X.7HN0  DATA  illX.IS) 
PTAB= (FLOAT < TLOW ) /SUM )*1  0  0.0 

CTAB  =  PTAF 

CUMMd  JaCTAS ~ 

WRITE  (6.12605  TLOW .PTAB* CTAB 
1260   FORMAT  ( 1 4X«9HBEL0W  LOW. 1 1 X. I 8«2 (3X.F1 2.4 ) > 


DO  350  1=1 .CON 

VALUE  =  DRANG(T> 

VAL1=DRANG( 1+1 ) 
PTAB=(TABE(  I  )/S'UM)*l    ,0 
CTAB=CTAB+PTAB 
CUMHf T+1)*CTAB 


350    WRITE  (6«1270)  VALUE«VAL1 »TA8E( U «PTA8»CTA8 

1270   FORMAT  UX«P12.4«6H   TO   «F12.4«3X«Fa«0»2<3X»F12.4> > 

I=C0N+1 

PTAB=(TAnF(  IjVSUM)*! 70 ~ ~~ ~~ '" 

CTAB=PTAB+CTAR 
CUMH(C0N+2)=CTAR 

WRITE  (6.1280)  VAL1 . TA3E ( I ) .PTAB « CTAB . SUM . CTAR . CTAB 

128U   FORMAT  UX«Ft2.4.l0H   OP   OVER  i  1  .X*F8.d.2Hx.Ft2«4) »/»  iSX«SHfOTAL 
1«  HX«F8»0»2<3X«F12»4n 

c        **#**»**********♦♦**»##**#*♦»##»*♦##♦♦## 

C    IFr..Il^yF>Sq.'»  1  |  PERFORM  HYPSOMETRIC  ANALYSIS _ 

T     ***»»»*»******»*»**»»»»»»»»»»»»»»»»»»»»» *~ 

IF  (IHYPSO)  450.451.450 
450    CALL  HYPSO(DRANG.ZMAX.ZMINiCON»CUMH,TlTL> 
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___    ______ 

END 
SIBFTC  STAT 

SUBROUTINE    gTA___________  I  I  .  JU.  X.ynAR  .YVAR  .  y<?  I  G  i  yMAY.yM  I  KM 

C  **********************  *****************  ************  ^**  ****--.>*********--,  ********** 

C       PROGRAM  TO  FIND  STATISTICAL  MEASURFS 

£  *  *  *  *  *  *  -if  *  *  *  *  *  *  *  *  *  *  *  *  *  -if-  *  *  -if-  *  -X-  *  *  *  -X-  *  *  *  *  *  *  *  -X-  *  -X  *  *  -"-  *  *  *  *  *  *  *  *  *  *  -"-  *  *  ■>*-  *  "-  - '-  *  *  *  *  *  *  *  * J '-  *  *  *  *  *  *  *  *  *  * 

INTEGER  QQV'F.COL? 

DIMENSION  y< 1  I « JJ) 

ENO"ROWS*COLS 

XSOP=0.0 
XSUM-0«0 

xmax=x( l . l ) 
xmIn=x(  l  .  l ) 

DO  10  1  =  1  »RO'"S 

DO  10  J=l  .COLS 

IF  <XC I . J) .GT.XMAX)  XMAX=XCI»J) 
IE  (X< I . J) .LT.XMIN)  XMIN=X<I«J) 
XSUM=XSUM+y( I • J) 

10     XSQR=XSQR+X< I .J)*X( I . J) 

XRAR=X?UM/FMO 

XVAR-  (xS0P-2i0*yDAR*XSU'«FMO*XDAP*xn/>P)/FMO 

X?IG=SOPT<XVAR) 

RETURN 


END 
SIBFTC  HYPS 

SUBROUT I NE  HYPSO ( DRANG i 7MAX . ZM I N ♦ CON  .  CUMH .  T 1  TL  ) 
C  ***-**»»*»********»**»***»****»»*»**»***  ^-^^^y^*.^.**.**.),.*^  **4<  .tJ.->.»»<t<.»»»»JH(.»#y^M(. 
~C     A  PROGRAM  TP  PERFORM  HYPSOMFTPIC  ANALYSIS* 
C ****** *********  **************  ****************************##->#*******_.(!.*.£.((..,(..)<..,..*;,.* 
C         REEERENCE- 

C Ai  N.  STRAHLER-195?-HYPSQMrTPIC  (AP^fi-ALTlTU^I  Ar.'ALY?  I  ". 

~~C  OF  EROSIONAL  TOPOGRAPHY-GEOL  •  SCC.  Amfq  1 C  A  BULL  •  VOL.*"? 

C  PAGES  1117-114?  .' 

C ************************************************** ************#***#*x**4»#*#**## 
DIMENSION  I  ROW (49. 81 ) «RELH(25) .RELA(?E) 

DIMENSION    TITL(P)  rctfMPT_7*l  . nPANG7?7> 

DATA  nOT.=LANK.cTAP/lHt | 1H  ilH*/ 

R^AL  IPOW 

Integer  con 

~c~       ***********  *****************  *********  ******** 

C    MAKE  PLOT  ALL  BLANKS* DRAW  BORDERS  WITH  DOTS .AND  PLOT  A  STAR  AT  TOP 
C    LEFT  AND  BOTTOM  RIGHT  CORNERS 

C     *******************************#************#  

DO  2  m=1 .40 


DO  2  N=l i 81 
IROW(MiN) =BLANK 

j_0_<l*l)=?TAR 


DO  3  N»2«81 
1R0WU «N)sDOT 
I  ROW (49*81 )=STAR 
DO  4  N-l .80 


I  ROW ( 49. N) -DOT 
DO  5  M-2.4R 
I ROW (M» 1 )=DOT 
I  ROW (M. 81  )=DOT 


C     **********************#******#************•*•*•♦ 
C    DETERMINE  POSITIONS  OF  STARS  WITHIN  THE  PLOT 
C     **************************##***##********♦#*♦ 
___N__ 


DO  20  1=1 .K 

RrLH( I )» ( DRANG ( I ) -7M IN ) / ( ZM AX-_M IN ) 

RELA  (  I  )  =  1  •  0—  C  (CUMH(  I  )-CUMH(  1  )  >/<CUMH (CON+1  )-«-llMH(  1  )  ) 

"  <49.0-<RELH(  I  )»49.  L__L_ 


N=(RELA( I >*81.0>+0«_ 
20      IROW(M.N)=STAR 

c        ****#***********##*******#♦#*#*♦♦*#******#*** 

C    PRINT  OUT  ALL  THE  PLOT  SEQUENTIALLY 


C     ****************************  *#******#.)|.******* 
WRITE (6*1 000)  TITL 

1000  FQPMAT( 1H1 .22HHYPS0METRIC  CURVF  FOP  .12A6) 
WRITE(6.1001  ) 

1001  FORMATdH  «8HREL»  HT«/) 
WRIT_(6«1002) 

1002  EORMATdH  .8X«  3H0.0.37X.  3H0  •5.37X.  3H1  .0  ) 
WRITE  (6.1003) 

1003  FORMAT ( 1 H  »  9X . 1 H+  «  39y . 1 H+ . 39y »  1  H+  ) 
WRITE(6. 1 004)   ( IROW< 1 »N) .N=l .81 ) 

1004  FORMATUH  .4X.BH1.0  +.81AI|5H+  1.0) 


48 

DO  100  M«2«24 


100  WR ITE < <M  1 005 )  <  tR0»i/<M»N)  «N=1  t81  ) 
1005   FORMATdH  «9X»P1A1)  " 

lr  n6   FORMATS  1H  (4y,RHn.=i  +|R1  A]  |BH+  aTWi " 

DO  101  M=26«4R 

101  VRITE(6i 10051  <IROW(M»N ) iNa 1 »81 ) 
WRITE (6< 1 007)  (IR0W(49»N)«N=ltflM 

1UU/   W^I11^00^,5HO•6  +'^A1^H+  ol6> 

WRITE<6«1002> 

WR1TE(6> 10081 

ICDS   FORMAT  (1H  ,44X«  1  3HRH.ATI  Vc    ARFA  /) " 

C       ***************************#*******.».#*.*.*.***** 

g  TH^^IT-Q^^N^^J^PA^r    USING    °^CIMAL    ViLU"    -    ^IMTS    POP 


AREA=0»0 

DO  200  I»l iCON 

^ A^A  =  AbFAlrRELH<  ™  '  >*<R^A<  '  >-^L.A(t  +  )  )  >  )/?.* 

WRITE <6i 1 010)  APFA 
1C10   FORMATdH  ,77HTHE  HYPSOMETRIC  INTEGRAL  fop  THp  Annvr  rn^/P  i     rue  a 
1REA  BELOW  THE  CURVF  )  t^  «F6i4)  "HOVE  CUR\/F  (  thc;  / 

~C      **********************»»»»4#»»i4»*».)i.».fc.it^ _ 

C    WRITE  OUT  COORDINATES  OF  HYPSOMETRIC  CUPVF 
C     *******************♦*»*#***#*#********♦* 

WRITE    (6*2001) 

guui       FOHWATdHl  «5ftHCALCULATFn    vALUFS    OF    HYPSOMFTPir    CURV    rnnDrMMATcr 

1«'^'}2^P0I^T    NUMBER.SX.SHA-AXIsTlOxTsHH-AXIS    ?                COORDINATES, 
g.1-...       yptTE     (6»2000)      ItRFLACI  )  iRfihI  I  ) 
^V™       g0RMAT<lH0,16X»!g.l  0X^6,44  10X>6. 4) — 

END 
SDATA 
20161611010,0  1 


M6F4?0?  MAP  °F  DRAINAGC;  5ASlN  »N  VERnuGO  H I LLS .  STPAHLEP  FIG  15  , 

9Q099099999Q9Q9999999oo99Q99PoooaooQQOQq  ,  -,  ,  =  J  ™ A-1  -T^e  I  a aSoSoSSo 
oqqoQoQaoooncooo^^™  "I"  Si?  7  =  9  ???  1  ^7*  1  fi!>R  61"  1  7  1  9  1  7-Qnnno 


llliigisiiililisiiiiiiiii- 

999Q,99Q9  I  31  21  230  U  551  166  1  ?  =  0  1  "?oonn,on0onQoonnnooconoa000-  tZZ22 

2SSSS222A  g§9L170!  l?1^!  ?qp 127^1  -^^Qooooooooooooooaooooooooo,-^  ,3,™ 

ll^^^U7511?01100117599PO9O9999O999aggQg95o^g^oo^|^||^°^ 


FIG.  C-3.  Typical  output  of   Program   GCON. 
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CONTOUR  MAP  OF  DRAINAGE.  BASIN  IN  VEROUGO  HILLS.  STRAHLER  FIG  IS  . 

320  OBSERVATIONS 

CONTOUR  LEVELS  SPECIFIED 

GRIDDED  DATA 

LINEAR  INTERPOLATION  WITMIN  GRID  SQUARES 

OBSERVED  VALUES 

C16F4.0) 

ROWS  =   20    COLS  =   16 

ZBAR  =    6145.009    7SIG  =    4262.532    MAxZ  «    9999.000    MlNZ 


1080.000 


VALUES  EMPLOYED 

ZMIN  »    1080.000    ZMAX  =    1880.000    CI 
GRID  *      *.    XINC  =       0.167    YTNC  = 
CON  a      16    PS  =       !3    COLMAX  a     91 
CONTOUR  SYMBOL  TABLE  IS  AS  FOLLOWS 

(ONLY  EVEN  SYMBOLS  PRINTFD  IN  BODY  OF  MAP) 


50.000 
0,278 


SYMwOL 


VALUE 


• 

NO 

DATA 

. 

BELOW   LOW 

0       1080 

.0000 

TO      1130. 

,0000 

1       1130 

.00(10 

TO      1180, 

,0000 

2      1180 

,0000 

TO      1230, 

,0000 

3      1230 

.0000 

TO      1280. 

,0000 

4        12HD 

.0000 

TO      1330, 

,0000 

5      1330 

.0000 

TO      1380, 

,0000 

6      1380 

.0000 

TO      1430, 

,0000 

7      1430 

.0000 

TO      1480, 

,0000 

8      1480 

.0000 

TO      1530, 

,0000 

9      1530 

.0000 

TO      1580, 

,0000 

A      1580 

.0000 

TO      1630, 

,0000 

8      1630 

,0000 

TO      1680, 

,0000 

C      1680 

,0000 

TO      1730, 

,0000 

0      1730 

.0000 

TO      1780, 

,0000 

E      1780 

.0000 

TO      1830, 

,0000 

F      1830 

.0000 

TO      1880, 

,0000 

•      1880 

.0000 

OR   OVER 

FREQUENCY  TABLE 

INTERVAL 

UNITS 

PERCENTAGE 

CUMULATIVE 

NO  OATA 

3081 

BELOW 

LOW 

0 

0. 

0. 

_ 1080.0000 

TO 

1130. 0000 

39. 

1.7751 

1.7751 

1130.0000 

TO 

1180,0000 

98. 

4.4606 

6.2358 

1180.0000 

TO 

1230.0000 

171. 

7.7833 

14.0191 

1230.0000 

TO 

1280.0000 

184. 

8.3751 

22.3942 

1280.0000 

TO 

1330.0000 

223. 

10.1502 

32.5444 

1330.0000 

TO 

13BO.0000 

230. 

10.4688 

43.0132 

1380.0000 

TO 

1430.0000 

230. 

10.4688 

53.4820 

1430.0000 

TO 

1480.0000 

213. 

9.6950 

63.1771 

1480.0000 

Tn 

1530.0000 

182. 

8.2840 

71.4611 

_  1530.0000 

TO 

1580.0000 

146. 

6.6454 

78.1065 

1580.0000 

TO 

1630.0000 

136. 

6.1903 

84.2968 

1630.0000 

TO 

1680. 0000 

88. 

4,0055 

88.3022 

__  1680.0000 

TO 

1730.0000 

84. 

3.8234 

92.1256 

1730.0000 

TO 

1780.0000 

67. 

3.0496 

95.1752 

1780.0000 

TO 

1830.0000 

34. 

1.5476 

96.7228 

1830.0000 

TO 

1880.0000 

17. 

0.7738 

97.4966 

1880.0000 

OR   OVFR 

55. 

2.5034 

100.0000 

TOTAL 

2197. 

100.0000 

100.0000 

FIG.C-3.  contd. 
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CONTOUR  NAP  OF  DRAINAOE  BASIN  IN  VEPDUOO  HILLS, STRAHLER  FIG  15 


>  e 

>  EEC 

>                                                                 EEEEE 

>                                                                  EEEEE 

1                                                                  EEEECE 

>                                                                  EEEEEC 

>                                                      CCCCCCCCC     EEEEE 

1                                                      CCCCCCCCCC        E 

CCCCCCCC 

>                                      aaaaaaaaaoaaaaaaaaaaa           CCCC 

AAAA                AAAAAAA                CCCC 

BH         AAAAAAA            CCCCC 

aaeeesas     aaaaaaaaa    CCCCCC 

-* 

HB8HH«flH8B«8e8             AAAAAA      CCCCCCCC 

fl8fl8BSR8eeeasfl8RR   aaBHaBBaasBB   aaaa    ccccccccccccc* 

aaaaaa  aaaaaseaaaaaae         anaasaaBB    a«aa        ccccccc 

>                 BaaRaaaaaaBaasea               Bsaaas   aaaaaaa 

ccc 

60B666666t>         888      AAAAAAAAAAAAAAAA 

6666666666666666      BBa       AAAAAA    AAAAAA 

66           666666666666        666    888                    AAAA 

6666666  66666666666666         6666   8B88B 

666666666666666666666            6666   88668888888888888 

66666666666666                A     666    88868888       88888 

6666666666                     ««     666     88             888888 

66666666                       AAA     6666                  888888 

66666        AAAAAAAAAAAAAAAAAAA    6666666666666666     8888888 

6           AAAAAAAAAAAAAAAAAAAAA           666666666 

6        AAAAAAAAA            AAAAAAAA4A        6666 

AAAAAAAA                AAAAAAAAA       6666 

• 

AAAAAA4A                   AAAAAAAA       6666 

666      AAAAA           22222222        AAAAA      66666 

6       AAAAA         2222222222       AAAAA       66E 

AAAAAA        22222222222       AAAAA      6666 

AAAAA       22222222222222     AAAA      6666666 

AAAAA       22222222222222     AAAA       6E 

AAAAA      222222222222222     AAAA       66 

•j 

AAAAA      7222222  2222222     AAA    6666666 
AAAAA     22222222   222222     AAC 

AAAA     22222       22222     AAA 

»■ 

AAA     2222           222     AAAA 

AA    2222            2222     AAAA 

• 

A    2222           2222     AAAAA 

•i 

222       0     2222     AAAAAA 

222        0     2222     C 

* 

222        no    22222 

•i 

222    0000000    2222 

222    0000000     22A 

22     00000000    22 

*« 

0000008000 

DOS 

0                                                                        • 

HYPSOHETRIC  CURVE  FOR   CONTOUR  NAP  OF  ORAINARE  BASIN  IN  VEROUOO  HILLSiSTRAHLER  FIS  IS 
ML.  HT. 

0.0  0,5 


0.0  0.5 

RELATIVE  AREA 
THE  HYPSOMETRIC  INTEnNAL  FOR  THE  ABOVE  CURVE  I  THE  AREA  RELOK  THE  CURVE  )  IS  0.4285 


